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Perturbative Wilson loops with massive sea quarks on the lattice 
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We present 0(g 4 ) calculations of both planar and non-planar Wilson loops for various actions 
in the presence of sea quarks. In particular, the plaquette, the static potential and the static self 
energy are calculated to this order for massive Wilson, Sheikholeslami-Wohlert and Kogut-Susskind 
fermions, including the mass and nf dependence. The results can be used to obtain ctjj-g and 
mb(rrTb) from lattice simulations. We compare our perturbative calculations to simulation data of 
the static potential and report excellent qualitative agreement with boosted perturbation theory 
predictions for distances r < 1 GeV -1 . We are also able to resolve differences in the running of the 
coupling between rif = 2 and n/ = static potentials. We compute perturbative estimates of the 
"/3-shifts" of QCD with sea quarks, relative to the quenched theory, which we find to agree within 
10 % with non-perturbative simulations. This is done by matching the respective static potentials 
at large distances. The prospects of determining the QCD running coupling from low energy hadron 
phenomenology in the near future are assessed. We obtain the result A^Lro = 0.69(15) for the two 
flavour QCD A-parameter from presently available lattice data where r^ 1 ~ 400 MeV and estimate 
a%{m z ) = 0.1133(59). 
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I. INTRODUCTION 

The calculation of Wilson loops in lattice perturba- 
tion theory is useful in a number of ways. An im- 
portant application is the prediction of a strong cou- 
pling constant a^L(mz) from low energy hadronic phe- 
nomenology by means of non-perturbative lattice sim- 
ulations ||, H %■> H Hi- Perturbative calculations of 
Wilson loops are also employed in the context of mean 
field improvement programmes of the lattice action and 
operators ||. 

In the limit of infinite Euclidean time separation, T — > 
oo, Wilson loops give access to the perturbative lattice 
potential. Our results improve the understanding of the 
violations of rotational symmetry at short distances of 
this quantity, as well as of the short distance effects of 
including various flavours of sea quarks. Furthermore, 
in the limit of large distances, R — > oo, the self-energy 
of static sources can be obtained from the potential, en- 
abling the calculation of rnb{rnt) from non-perturbative 
simulations of heavy-light mesons in the static limit || . 

Perturbative computations of Wilson loops on the lat- 
tice now reach back as far as two decades. For this rea- 
son, much of the groundwork for the calculation covered 
in this article is well known, and has been reproduced by 
us. In addition to improving the precision of many pre- 
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viously existing results and presenting them in a more 
comprehensive way, we extend the calculations by incor- 
porating the actions, sea quark masses and observables 
that are relevant for present day lattice simulations. We 
also perform the first calculation of non-planar Wilson 
loops. Because of the degree of replication we present a 
brief survey of the existing literature here. 

In 1981 Miiller and Riihl |1(| performed a highly rig- 
orous calculation of the 0{<p) contribution to the Wil- 
son loop for pure U(l) and SU{2) gauge theories and 
presented their result in terms of lattice integrals. In 
particular, they obtained a closed form for the T — > oo 
limit and hence the static potential. The plaquette in 
SU(N) pure gauge theory was first calculated to this or- 
der by Di Giacomo and Rossi ]lT| l , a gain in 1981. In 
the same year Hattori and Kawai jl^] estimated 0(g 4 ) 
contributions to larger Wilson loops, by modelling some 
of the relevant lattice contributions by their continuum 
counterparts. Rigorous results to this order were subse- 
quently independently obtained by three groups of au- 
thors: Weisz, Wohlert and Wetzel Q, Curci, Paffuti 
and Tripiccione and Heller and Karsch jl5|. The 
first reference is the most general and also applies to 
Symanzik-improved pure Yang-Mills actions. However, 
numerical results on Wilson loops are only available for 
the Wilson action. The 0(g 4 ) fermionic contribution to 
the plaquette was obtained by Hamber and Wu jl6) for 
Wilson quarks and for Kogut-Susskind (KS) quarks by 
Heller and Karsch p7t| . More recently, Panagopoulos et 
al. calculated the 0(jf) SU(N) gluonic contribution Q 
as well as the 0(g 4 ) and 0(g 6 ) corrections [|l9| in pres- 
ence of massive Wilson fermions. The only numerical 
calculation to-date of small Wilson loops for Symanzik 



improved gluonic actions was performed by Iso and Sakai 
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More recently, the Heller and Karsch calculations |15[| 
have been extended by Martinelli and Sachrajda || in 
a manner similar to that of Miiller and Riihl flCf] : the 
R, T — > oo limit of the perimeter term of Wilson loops 
has been taken to obtain the static self energy for SU (3) 
gauge theory, with the addition of massless Wilson and 
Sheikholeslami-Wohlert (SW) sea quarks. 

The structure of this paper is as follows. In Sec. || we 
describe the procedure used to calculate Wilson loops in 
perturbation theory and present the perturbative expan- 
sions of small Wilson loops, both planar and non-planar, 
to 0(g 4 ) using a variety of actions. Although partially 
addressed by previous authors jl3[ [l^, the static 
potential has not been well treated in literature, which 
we attempt to rectify in Sec. where we also obtain 
the static source self energy from the asymptotic R — > oo 
potential. This sets the stage for Sec. IV where these 
perturbative results are confronted with state-of-the-art 
lattice data from simulations with Wilson, SW and KS 
fermions: we determine the "/?-shift" and the running 
coupling ajjg, before comparing perturbative and non- 
perturbative potentials at short distances. In Sec. |v| we 
conclude and summarize. Our article is augmented by 
two Appendices on the relations between different lat- 
tice and continuum renormalization schemes and on the 
perturbative /3-shift. 




FIG. 1: Tree-level [0{g 2 )] contribution to the Wilson loop. 
The coordinate space propagator is convoluted around the 
possible source and sink links of the loop. This is done by 
constructing the corresponding sum of phase factors in mo- 
mentum space and Fourier transforming. 



extent of T lattice units pointing into the time (4) direc- 
tion, that are separated by a spatial distance Ra. The 
smallest such Wilson loop is an a x a square, the plaque- 
tte □ = W(l, 1) which, appropriately normalized, is the 
expectation value of the Wilson gauge action, 



-P E -^R-etrU^, 



(2) 



II. WILSON LOOPS IN LATTICE 
PERTURBATION THEORY 

In this section we explain our method and notations, 
before displaying results on small Wilson loops with mas- 
sive fermions to 0(g ). We also include some findings 
that apply to the Iwasaki and Symanzik-Weisz gauge ac- 
tions. 



where (3 = 2N/g 2 and U x ^ v = U x ^U x+a p, tV Ul +aj> ^Ul >u . 
Unless stated otherwise this is the gauge action which we 
combine either with the Wilson, SW or the KS fermionic 
actions. 

Writing Ui = e l9aAl , where I = (x, /i) denotes the link 
connecting x with x + aji and Ai is a short hand notation 
for An(x + §/i), and expanding the exponentials within 
Eq. (|l|) one obtains, 



W{C) = 1 - g 2 W 2 - g 3 W 3 - g 4 W A 



(3) 



The Method 
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We consider SU(N) lattice gauge theory on an infinite 
four dimensional isotropic hyper-cubic lattice with lat- 
tice spacing a. An oriented closed curve C of length na 
touches a sequence of sites, {xi : Xi/a S Z 4 } such that 
Xi + i —Xi 6 {±afi : )jl = 1, . . . , 4}, x n +i = X\. The Wilson 
loop W(C) around such a curve is the expectation value 
of the path ordered product of gauge links, 



n u t 
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(i) 



where /i^ G {±1, 
by x i+ i - x. 

We write W(R, T) for a "rectangular" Wilson loop 
where the curve C contains two opposite lines with an 
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W2 = ^ TaTb \L A t A i 

\h,h i 



W 3 



±4} denotes the direction indicated W& 

t 



6N 



tr T a T b T c I A h A h A i : 



3 £ [AlAlAl-AlAlAl] 



24N 



E 

Jl<i2<'3 



li<h<h 



tTT a T b T c T d I 24 



E 



(4) 



(5) 



Af A) Af Af 

ll l>2 t3 U 



\ /l</ 2 <^3<U 

(AlA b h AlAi 3 +AlAlAlAi 



3 



+ AlA\ 2 Ai 3 Ai) 



(6) 



4A? 1 4 2 ^ 2 4 2 ; 



Y^AlAlA^Ai 



k 


1 2 

ll 


k pj 

: P >sW 

k 


k 


1 2 

ll 


\ / 1 

k 


li 




b 



FIG. 2: Pure gauge one-loop [C(g 4 )] vacuum polarization 
contribution to the Wilson loop. This is the sum of continuum 
like three-gluon and ghost-ghost-gluon couplings, and both 
a continuum like four-gluon tadpole contribution and lattice 
specific four-gluon and two-ghost-two-gluon tadpole pieces. 
The last measure diagram is due to the Jacobian involved in 
changing variables in the path integral from Ui to Ai. 





FIG. 3: The leading rif contribution arises at 0(g 4 ) through 
the fermionic contribution to the vacuum polarization, and 
the lattice specific fermionic tadpole contribution. 

The expectation values in the above formulae depend 
on g 2 and are sums of gluonic position space n-point 
functions. The leading order [0(g 2 )] contribution to the 
Wilson loop is the sum of tree-level two-point functions 
depicted in Fig. [I]. This is easily expressed in momentum 




FIG. 4: The 0(g 4 ) spider contribution to the Wilson loop. 
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FIG. 5: One-loop [0(gr 4 )] two-gluon exchange contribution 
to the Wilson loop. One part of this diagram is accompa- 
nied by a Cp/2 colour factor, that arises when exponentiating 
the tree-level contribution. The remainder is proportional to 
CfN/2 and specific to non-Abelian gauge theories. 



space as a sum of phase factors multiplying the momen- 
tum space propagator: 



W 2 
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For the Wilson gauge action the gluon propagator reads, 
S ab (S - S ) 



(^(-fc)4 2 (fc)) 



b fu\\ — u vW>P2 



k 2 



IT 



ah 



Xk), 



(8) 

where k ^ = 2a _1 sin (k^a/2) and II is of order g 2 . The 
colour trace then results in the Casimir factor of the 
fundamental representation, Cf = jftr T a T a = (TV 2 - 
1)/(27V). All formulae in this Section as well as in Sec. [II 
below apply to the general SU (N) case unless stated oth- 
erwise. In addition the constants Cf are factorized out 
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in such a way that the results apply to external sources 
in any representation D of the gauge group, under the 
replacement Cf * Co- 

In explicitly summing the above phase factors one eas- 
ily obtains the well known tree-level expressions for pla- 
nar Wilson loops. The vacuum polarization tensor II has 
to be considered to leading order for an 0(g 4 ) evalua- 

(2) 

tion of Wilson loops which requires W 2 ■ The relevant 
graphs for the gluonic sector are depicted in Fig. 0. The 
fourth diagram involving a two-gluon-two-ghost vertex 
(which is lattice specific) as well as a lattice contribu- 
tion to the four-gluon vertex of the third diagram, result 

(2) 

in a lattice tadpole correction to W 2 that is exactly 
(2N 2 - 3)/(24iV) times the tree-level value W [ 2 Q) . The 
two additional diagrams of Fig. [| have to be considered 
in the case of dynamical fcrmions. 

To the order at which we are working throughout this 
paper, W3 only contributes through its contraction with 
the three gluon vertex, depicted in Fig. |], and we write 



W 3 = gW. 



(i) 



In momentum space the phase factor associated with a 
typical entry in the sum Eq. (||) is 



^_2^_2__ irT a T b T c ik[x 2 -x!+^ (A2-A1)] 

xe ife'[x 3 -x 1+ § (A3-A1)] (A a ^{-k - k')A b (k)A° 3 (&')}■ 



The term 



W± = wi 0) +0{g 



(11) 



0(g 3 ). 



(9) 



only contributes to the Wilson loop at 0(g 4 ) through 
its contraction with double gluon exchanges, as depicted 
111 Fig. § The generic phase factor associated with this 
diagram is given below, where due to the lack of an ap- 
propriately symmetrized contracting vertex we perform 
the Wick contraction manually, preserving the ordering 
of colour factors: 



" /a d 4 k d 4 k' i 
_„ /a (2tt) 4 (2tt) 4 TV 



2 ifc[x 2 -xi+f (£2-/4i)] e ifc'[x4-x 3 +§ (fii- fo)) lj±a (—jfc)^ 6 (J^j\c 
ife[x 3 -x 2 + f (A3 -p,2 )] gifc' [x4 -xi+f (A* - Ai )] [—k')A h (—k)A c ^ (fc)yl^ (fc')) 



(12) 



After evaluating colour traces, this term can be split into 
an "Abelian" and a "non-Abelian" part. The Abelian 
part is proportional to Cp/2 and can be obtained by ex- 
panding exp(— g 2 W 2 ). The remainder is specific to non- 
Abelian gauge theories and proportional to CfN/2, the 
same colour factor that also accompanies the vacuum po- 
larization diagrams. 

The integrands were coded with the help of a C++ 
package written by one of the authors, and then inte- 
grated using VEGAS p2|] . For finite Wilson loops to 
0(g 2 ) and 0(g A ) we compute 4 and 8 dimensional mo- 
mentum integrals, respectively. In the case of the poten- 
tial one can analytically perform the limit T — > 00. This 
allows us to compute the momentum integrals in one less 
dimension, leaving 3 and 7 dimensional integrals. 

The fermionic Wilson-SW and KS vacuum polariza- 
tion graphs have been computed after a further analyti- 
cal integration over the 4-component of the internal loop 
momentum, leaving a 6 dimensional integration to be 
performed in the cases of the potential and the static self 
energy. 

It is relatively easy to automate the generation of the 
phase factors that are required for arbitrary curves C. 
We generate the perturbative expansion for both small 
rectangular Wilson loops and for the non-planar loops 
shown in Fig. 0. We refer to these non-planar loops as 



the chair and the (three-dimensional) "parallelogram" ; 
both are contained within the unit cube and are used 
within certain improved gauge actions. 




FIG. 6: The non-planar Wilson loops in the unit cube which 
we refer to as the chair and the "parallelogram" , respectively. 



B. Definitions 

We write the tree-level contribution to the expansion 
of the Wilson loop with the Wilson gauge action as, 

M/ 2 (0) = C F W T - (13) 
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For Wilson and SW fermions with the Wilson gauge ac- 
tion we write the vacuum polarization insertions as, 



W, 



(2) 



'N 
Cr ( y Wn 



27V 2 



24A 



-Wt + nfXf 



For Wilson- SW fermions we define, 

A / — A / + A / c sw+^f c sw , 



(14) 



(15) 



where c$w = in the Wilson case and csw = 1 + C(ff 2 ) 
in the SW case. The Xf will depend on the quark mass 
in lattice units ma. In the case of Wilson-SW fermions 
and the order of perturbation theory we are working at 
we have, ma = moa = (ft -1 — k~ 1 )/2. As anticipated 
above, we have factorized out a lattice tadpole term that 
is proportional to Wt- For KS fermions we define anal- 
ogously, 



X 



= X { f KS) . 



(16) 



In this case Uf has to be a multiple of four. For each of 
the above quantities we may write a similar expansion in 
terms of «l = <7 2 /(47r), rather than g 2 , with lower case 
coefficients w and x, where wt = iirWr, wn = (47r) 2 Wn, 
x = (47r) 2 X and so forth. 

We denote the spider contribution by, 



W. 



(i) 



C F N 



W s . 



(17) 



The double gluon exchange can be written as the sum of 
Abelian and non-Abelian pieces, where the Abelian piece 
is the next term in the exponentiation of the tree-level 
contribution, 



colour sources, D, through its proportionality to the cor- 
responding Casimir charge Cjj — (tr TgTp) / 'diron . Note 
that in the case of the fundamental representation, by 
using the relation C F = (N 2 - 1)/(2N), Eq. (§]]) can be 
rearranged into the form, 



1 



where 



W NhQ = (N z - 1) ( W a + —W b 



W a = — (6W S + W T ~ 3W T ) , 



W h 



1 

16 



-W T + 2W T ) 



(23) 

(24) 
(25) 



For improved gluonic actions Eq. ( |2l| ) does not apply 
and hence W a and Wb can take somewhat different 
forms. The factorization Eq. ( p3| ) is frcouently employed 
throughout the literature, e.g. in Refs. |1^, |l^, [2p[| . 



C. Small Wilson loop results 

The pure gauge results on small Wilson loops for the 

Wilson gluonic action are displayed in Tab. |j. We re- 

produce_the known plaquette results of Di Giacomo and 

l 



Rossi Jri) and Alles et al. 1 19 1 and the small planar 
Wilson loop results of Wohlert, Weisz and Wetzel [fT3| . 
Boldface values have been calculated by us for the first 
time. We also calculated tree-level Wilson loops for the 
Symanzik- Weisz (SyW) Q and Iwasaki (I) [2(| actions. 
These results are displayed in Tabs. |l| and fflj respec- 
tively, together with the 0{g A ) values of Iso and Sakai |2^| 
(italicized) which we have not validated ourselves. 



(0) 



r 2 



C F N 



W NA - 



(18) 



In combining the above Eqs. ( |l3| ) - ( ^8[ ) with Eqs. (g), 
(fjj), (||) and (pd|), we obtain the expansion of the Wilson 
loop, 



W = l-W Jj0 g 2 -(W Nh0 + C F n f X f )g i , (19) 



where 



W LO = C F W T , 

( N 2N 2 - 

W^o = C F (-W* + _ 

Wy, = Wu + Wna + Ws 



*W T -^W T 



(20) 

(21) 
(22) 



The W T term within Eq. (|2l]) originates from the expo- 
nentiation of the Abelian part of the one gluon exchange, 
Eq. (p^|). It is cancelled in the Taylor expansion of InVF, 
which turns out to be proportional to C F , at least to 
0(g ). The latter observation, which is referred to as 
the "Casimir scaling" hypothesis in the literature (see 
e.g. Refs. [^3[ ||)), implies that, for given N and n^, 
InW only depends on the representation of the static 



TABLE I: Pure gauge Wilson loops with the Wilson gau ge 
action. The different contributions are defined in Eqs. (n9|) - 
(p2[). The values VKnlo a PPly to Wilson loops in the funda- 
mental representation of 5(7(3) gauge theory. Boldface values 



are calculated here for the first time. 


loop 


W T 


W s 


Wnlo 


1 x 1 


0.25 


0.0100109(4) 


0.033911(1) 


chair 


0.3922(1) 


0.02204(2) 


0.01629(1) 


par all. 


0.4267(1) 


0.02730(2) 


0.01128(1) 


1x2 


0.43110(6) 


0.02463(2) 


0.00382(1) 


2x2 


0.68466(8) 


0.05303(2) 


-0.12043(9) 



We are able to reproduce the known contributions to 
the plaquette for massless KS |l7| and Wilson Q quarks 
with increased precision. The calculation of the SW con- 
tribution 2 , as well as all results for massive quarks are 



1 In fact the result of this reference is more precise than ours: 
they obtain W^LO = 0.033910993(1) while we quote W NLO = 
0.033911(1) in the table. 

2 Our result has already been used in Refs. 
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TABLE II: Pure gauge Wilson loop results with the 
Symanzik-Weisz action. Italicized values have been obtained 
in Ref. W a and W b are defined in Eq. (§|) and Wnlo 

applies to TV = 3. 



loop 



W, 



SyW 



W SyV 



w, 



SyW 



SyW 



NT.O 



1 x 1 0.18313(1) 

chair 0.28850(8) 

parall. 0.3162(1) 

1 x 2 0.33130(6) 



-0.001133 0.001333 



-0.00678 



0.00830 



-0.00788 



-0.0469 



new. We also calculate the one- loop n/ piece of the 
plaquette with the Iwasaki gauge action which has been 
used, for instance, in the dynamical fermion simulations 
of the CP-PACS grou p |(| . We display a selection of 
vn == results in Tab. 



IV 



In Fig. [?] we show the mass 
dependence of Xf for the case of the Wilson gluonic ac- 
tion, combined with all three different fermionic actions. 
The open symbols correspond to the respective m = 
limits. The numbers that are plotted in the figure are 
also displayed in Tab. 0. 

TABLE III: The same as Tab. |n|ior the Iwasaki gauge action. 



loop 



Wl 



NLO 



1 x 1 0.10514(2) 

chair 0.16676(4) 

parall. 0.18494(4) 

1 x 2 0.20166(4) 



-0.002269 0.003142 -0.01536 



-0.00653 0.00881 -0.04441 



TABLE V: Mass dependence of the fermionic contribution to 
the plaquette for the Wilson gluonic action. 



ma 


x) 710 


X) 710 


x) 710 


x) 710 


0.00 


-1.2258(7) 


-1.392(3) 


0.0404(5) 


-1.1927(2) 


0.05 


-1.2211(7) 


-1.275(2) 


0.0267(5) 


-1.1769(2) 


0.10 


-1.2076(6) 


-1.176(2) 


0.0154(5) 


-1.1603(2) 


0.20 


-1.1576(6) 


-1.006(2) 


-0.0011(4) 


-1.1256(2) 


0.50 


-0.8954(5) 


-0.655(2) 


-0.0207(3) 


-1.0158(2) 


1.00 


-0.4617(4) 


-0.359(2) 


-0.0208(2) 


-0.8455(1) 


2.00 


-0.1091(2) 


-0.141(1) 


-0.0104(1) 


-0.5967(1) 


4.00 


-0.01212(6) 


-0.038(1) 


-0.00261(4) 


-0.33661(4) 


8.00 


-0.00093(2) 


-0.0070(6) 


-0.00034(1) 


-0.14920(2) 












t ■ » - 


-0.5 






-1 














□ 




-1.5 






-2 




KS ' - 






Wilson ■ 






SW 


-2.5 







0.03 0.1 0.3 



1 3 10 30 
ma 



TABLE IV: Massless fermionic contribution to Wilson loops 
for the different quark actions. The plaquette is displayed for 
both Wilson (W) and Iwasaki (I) gluonic actions. All other 
loops are just shown for the Wilson gauge action. Xf are 
defined in Eqs. (^ and (|li|. 



loop 


X ( f KS) /W~ 3 


x< 0) /io- 3 


X^/IO" 3 


AfVlO" 3 


W: 1 X 1 


-1.2258(7) 


-1.392(3) 


0.0404(5) 


-1.1927(2) 


I: 1 x 1 


-0.2592(4) 


-0.294(2) 


0.00156(2) 


-0.3396(2) 


chair 


-2.674(2) 


-2.974(7) 


0.0724(9) 


-1.9672(6) 


parall. 


-2.890(2) 


-3.31(1) 


0.101(1) 


-2.3248(6) 


1x2 


-2.357(1) 


-2.652(6) 


0.113(1) 


-2.9518(6) 


2x2 


-4.113(3) 


-4.76(2) 


0.322(2) 


-6.684(1) 



FIG. 7: n/ contributions to the one-loop plaquette for KS, 
Wilson, and SW quarks. All Xf interpolate smoothly between 
the massless values (open symbols) and zero as the quark mass 
is increased. 



mass dependence remains weak for any Wilson loop with 
linear dimensions that are small compared to the inverse 
quark mass. 



III. THE PERTURB ATIVE STATIC 
POTENTIAL 



In the region ma < 0.25, which is most relevant for lat- 
tice simulations, this dependence is relatively weak and 
can be parametrized by, 

Xf(ma) — Xf(0) + bma + c(ma) 2 , (26) 

with Xf (0) as in the first row of Tab. |y| and b w 0, 0.0024 
and 0.0024 and c w 0.0017, -0.0023 and -0.0017 for KS, 
Wilson and SW fermions, respectively. As expected the 



After a historical survey of the existing literature we 
will briefly explain our method and the notations that we 
adopt, before presenting results on the 0(g ) potential. 
We will then discuss limiting cases, including the contin- 
uum limit and lattice artefacts as well as mass dependent 
terms in the conversion between lattice and MS schemes 
at finite lattice spacing. We conclude with paragraphs 
on the perturbative /3-shift, boosted lattice perturbation 
theory and the static source self energy. 
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A. The method and definitions 

The one-loop potential has first been calculated for 
pure SU(2) gauge theory by Miiller and Riihl Jk}. Sub- 
sequently, a closed formula for pure SU (N) gauge theory 
with the standard Wilson as well as Symanzik improved 
gauge actions was derived by Weisz and Wohlert |fl3| . 
The coefficients of the perturbative expansion have first 
been evaluated numerically by Altevogt and Gutbrod Q 
for SU(2) pure gauge theory with Wilson glue on 
isotropic and anisotropic lattices. 

Here we restrict our calculations to the case of the 
Wilson gauge action on an isotropic lattice but include 
massive Wilson, SW and KS fermions as well as N > 2. 
We also improve the numerical precision with respect to 
earlier results and incorporate off-axis lattice separations 
R of the colour sources. We define the static potential, 



aV(Ra) 



lim 



d\nW{R,T) 



i x dT 

2 , ir /T>\„4 



Vi(R)g 
vx(R)a L 



V2(R)g 4 
-u 2 (R)o;|, 



(27) 

(28) 
(29) 



where a L = 2 /(4tt), ^(R) = (4tt) 4 ^(R). W{R,T) de- 
notes a (generalized) rectangular Wilson loop with the 
temporal extent of T lattice units. The expansion of 
the potential in terms of g 2 is suitable for calculations 
in lattice perturbation theory and for the comparison 
with data from lattice simulations, while the expansion 
in terms of qj, turns out to be more convenient to relate 
our results to those obtained in a continuum scheme. 

The only gluon exchanges that contribute to the T 
dependence of W(R, T) in the limit T — >• oo [at least up 
to 0(<? 6 )] are those between temporal lines of a Wilson 
loop. After exploiting translational invariance in time, 
Eqs. (0) and (|) result in, 



V 1 (R) = C F 



d\ 2sin 2 (^ 
(2tt) 3 q 2 



(30) 



The 0(g A ) coefficient V2 consists of two parts: one con- 
tribution from W2, i.e. the vacuum polarization tensor 
II, inserted into Eq. (f30|), and a second contribution 
that originates from the "non-Abelian" part Wna of W4 
[Eq. dig)]. The spider diagrams, W3, do not contribute 
to V(R) at this order. A closed formula for V2 can be 
found in Ref. [ fi3|| , and we extend this calculation by in- 
corporating the fermionic contributions to the vacuum 
polarization tensor. In our calculation we produce V2 au- 
tomatically from the Feynman rules and find agreement 
with the published analytic form. 

The potential can be factorized into an interaction part 
and a part that is due to the static source self energy, 
V s = 2<5ra s tat, 

V(Ra) = Vint(Ra) + Vs. (31) 
Since in perturbation theory Vint — > as r — > 00 we can 



TABLE VI: Gluonic contributions to the static potential, [see 
Eqs. ([34J) - (p5|)]. The last line (R = 00) contains the respec- 
tive contributions to the static self energy aVs- The values for 
V2 in the last column only apply to the case N = 3, n/ =0. 
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(i,l,uj 
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n oonc / a \ 
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(1,1,1) 
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1.5697 (2) 


("1 Ef A At} / r\\ 

0.5443 (2) 


1 n noTO / c \ 

10.0378 (5) 


(2,1,0) 


01 coc\c\ 
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1 TO CT O / 0\ 

1.7253 (3) 
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1 7007 /'OA 
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A 7*30 /) 0^ 
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11 OOOO /'o\ 
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/O O 0\ 
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11. /2y2(lU) 


(3,3,0) 


23.4016 


1.9153 (5) 


0.8776 (4) 


12.0863(12) 


(4,2,0) 


23.4889 


1.9266(52) 


0.8914(11) 


12.1607(81) 


(4,2,2) 


23.6538 


1.9530 (6) 


0.9186 (4) 


12.3137(15) 


(3,3,3) 


23.7512 


1.9679 (7) 


0.9345 (4) 


12.4023(15) 


(4,4,0) 


23.8686 


1.9884 (7) 


0.9546 (4) 


12.5162(16) 


(4,4,2) 


23.9512 


2.0011 (8) 


0.9695 (4) 


12.5944(17) 


(4,4,4) 


24.1286 


2.0310 (9) 


1.0018 (5) 


12.7675(20) 


OO 


25.2731 


J2 ■ 3.5459(15) 


14.1122 (31) 



identify the self energy Vs with lim r _ >0O V(r) and write 

aVs = V 1 (oo)g 2 + V 2 (oo)g 4 + ---, 
aM„t(Ra) = V L , 1 {R)g 2 + V L , 2 {R)g 4 + -- 



(32) 
(33) 



where Vt,i(R) = V^(R) - V<(oo). The limit i? 00 can 
easily be realized by replacing the factors sin 2 [qRa/2] = 
[1 — cos(qRa)]/2 within the external 3-dimcnsional 
Fourier transformations in the expressions for V\ [cf. 
Eq. (p?j|)] and V2 by the constant value 1/2. 
We write, 



V(K) 
V 2 (R) 



C F V T (R) 

r n 

C F l-[V n (R) 



V NA (R) 



(34) 
(35) 



2A 2 



24N 



-V T (R)+n f V f (R) 



where Vu denotes the gluonic vacuum polarization and 
Vna is a contribution, specific to non-Abelian gauge the- 
ories, where the ordering of gluon vertices along the Wil- 
son loop has to be considered. In the case of KS quarks 
we denote the fermionic vacuum polarization contribu- 



tion to the potential by V 



(KS) 



f 



while for Wilson- SW 
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TABLE VII: Fermionic contributions to the static potential 



V f [see Eqs. (|3f 
Wilson-SW [F f w 



) and for massless KS [V) 

see Eq. (p6b] quarks. 



(KS)-: 



and 



TABLE VIII: Fermionic contributions to the static potential 
for massive KS and Wilson-SW quarks. 



R 


v ( f KS) /w- 3 


v f (0) /io 


—3 


vfVio- 3 


v f (2) /io 


-3 


(1,0,0) 


-1.0830 (1) 


-1.1318 


(!) 


0.0563 (2) 


-1.5070 (1) 


(2,0,0) 


-1.4893 (3) 


-1.6698 


(1) 


0.1469 (3) 


-2.6859(16) 


(3,0,0) 


-1.7004 (5) 


-1.9272 


(2) 


0.1978 (3) 


-3.0033 


(3) 


(4,0,0) 


-1.8123 (7) 


-2.0723 


(2) 


0.2331(34) 


-3.1578 


(3) 


(5,0,0) 


-1.8900 (9) 


-2.1659 


(2) 


0.2479 (4) 


-3.2589 


(3) 


(6,0,0) 


-1.9427(11) 


-2.2328 


(3) 


0.2611 (4) 


-3.3294 


(4) 


(7,0,0) 


-1.9851(13) 


-2.2818 


(3) 


0.2706 (4) 


-3.3840 


(4) 


(8,0,0) 


-2.0155(16) 


-2.3207 


(3) 


0.2779 (5) 


-3.4234 


(4) 


(1,1,0) 
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(!) 
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(1,1,1) 
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(!) 
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(2,2,2) 
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(4) 


(o o n\ 
^o,o,U ) 
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U.zoDD \0) 
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(4) 
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-1.8499 (9) 


-2.1233(11) 


0.2389(10) 


-3.1917 
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-1.8824 (9) 
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(4) 


(3,3,3) 


-1.9005(10) 
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-3.2536 
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CO 


-2.3359 (4) 


-2.6808 


(4) 


0.3266(10) 


-3.7174(12) 



quarks we split this contribution into three parts, 
V f (R) = Vj a) (R) + c S wV f m (R) + c 2 sw vj 2) (R). (36) 



B. Results on the potential and violations of 
rotational symmetry 

In Tab. [yj we display Vr, Vjj and Vna as well as Vi 
for N — 3, rif — for small on- and off- axis distances 
R. At the origin all V are zero and Vr(l,0, 0) = 1/6. 
In Tab. VII we display the corresponding results on th e 
different V/'s for massless quarks while in Tabs. VIII - 
IX the results for distances R < 3 as well as for R — > oo 
are shown for quarks of various masses. 

In Fig. ||we separately display the different C(g 4 ) con- 
tributions to the quantity — i?V^(R) as a function of R. 
The factor —R results in a cancellation of the leading 
order Coulomb behaviour. Note that the lattice tadpole 
contribution to the vacuum polarization (solid circles) is 
numerically equally important as the sum of the remain- 
ing vacuum polarization diagrams and the non-Abelian 
contribution. As one would naively expect, the effect of 
fermions (open squares) goes into the opposite direction, 
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0.25 
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relative to the pure gauge contributions. In a continuum 
calculation using dimensional regularization — i?V^(R) 
contains only a logarithmic dependence on R, but in the 
lattice calculation the R dependence of the points is dom- 
inated by a linear piece arising from the static source self 
energy contribution to V2(R), Vs,2 — ^2(00). In Fig. ^ 
we subtract this term, before multiplying the result with 
— R. Note the logarithmic scale. The static source self- 
energy was previously known to order a 2 with massless 
Wilson and SW quarks J9| and to 0(a 3 ) in the pure gauge 
theory |3^] . The mass-dependence, as well as the re- 
sults for KS sea quarks, are new. 
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TABLE IX: Tab. VIII continued 
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FIG. 8: 0(g 4 ) contribution to the static potential, multiplied 
by —R for N = 3. The solid triangles are the sum of the pure 
gauge contributions: non-Abelian, vacuum polarization and 
the lattice tadpole contribution to the vacuum polarization. 
The open triangles incorporate 2 flavours of Wilson fermions. 





0.14 




0.12 




0.1 


DC 


0.08 


CNJ 




> 


0.06 


DC 






0.04 




0.02 









-0.02 



sum: pure gauge 

glue w/o latt. tad's 
lattice tadpoles 
2 Wilson quarks 




10 



FIG. 9: The same as Fig. ^ but with the self energy con- 
tribution subtracted, V2,l(R) = V^(7?) — V2(oo), and with a 
logarithmic scale. 



In the limit r 3> a, i.e. R = r /a 1 rotational sym- 
metry should be restored such that we are able to com- 
pare our result to calculations performed in a continuum 
scheme (continuous curves in Figs . p| a nd ^). The po- 
tential has been calculated in the MS scheme to order 
a 3 |33|, [34j. The conversion between the MS and the 
lattice scheme has been worked out to this order too for 
Wilson-SW fermions |35|, |36| (for KS fermions only to 
0(a 2 ) H?], [jH), and Vs vanishes by definition to all or- 
ders in dimensional regularization. The f orm of the large 
R expectation is worked out in Appendix A 4 [Eqs. ( A72| ) 
- (A74)]. Here we display the parametrization of the 



curves in the limit of massless sea quarks: 



rV Lt2 (r/a) 



r I a- 



-rV Cj2 (r) 

47T 



(37) 



[61+ of + 2/3 ln(r/o)] 



The constants 61 and a{ = ai + 2^(3$ are defined in 
Eqs. (A14) and (A34), respectively. b\ originates from 



the c onve rsion between the MS and the lattice scheme, 
Eq. ( A13| ). The logarithmic running is proportional to 
the /3-function coefficient /3 = (UN - 2n/)/(127r). 
We can define an effective Coulomb coupling CfOcr 
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from the potential: 

-r[V(r)-V s ] = CpanXr- 1 ) (38) 
= C F a L - RV L , 2 (R)g 4 . (39) 

Note that in this limit, Cfcul = ~RVl,i(R)- In Fig. ^ the 
enhancement of this effective coupling, CfCur, is demon- 
strated, relative to the tree-level value as a function of 
R. The points are indeed consistent with the logarith- 
mic running proportional to 2(3qCf In R /(4n ) that is ex- 
pected from the QCD /^-function, Eq. (Al) [as well as 
from Eq. (37) above]. The lattice tadpole terms do not 
contribute to this running but renormalize the overall 
value of the coupling. Lattice simulations are usually 
performed around g 2 1, where a fit to quenched lat- 
tice data H © Y ields e ~ C F a R {l/(3a)] w 0.3. The 
tree-level expectation in this case, however, is substan- 
tially smaller: Cf&r — Cf&l ~ 0.106. As one can read 
off from the figure, the one-loop value at R = 3 adds 
about 0.08 to this, but still the non-perturbative result 
is underestimated by perturbation theory in terms of the 
bare coupling g 2 by more than 30 %. We will discuss the 
improvement resulting from the u se of so-called boosted 
perturbation theory in Sec. IV C below. 




FIG. 10: Visualization of lattice artefacts. Vl,i and Vl,2 
denote the tree-level and 0(g 4 ) lattice potentials, where the 
self energy pieces V;(oo) have been subtracted. V c ,i are the 
respective large R expectations. 

At small R values the lattice results scatter around 
the continuous curves. We shall investigate these viola- 
tions of rotational symmetry in more detail: in the top 
half of Fig. 10 we plot ratios of lattice and continuum 
C((/ 2 ) (open squares, Vl,i/V c ,i — 1) and 0(g 4 ) (open cir- 
cles, Vl,2/V c ,2 — 1) contributions for on-axis as well as for 
some off-axis distances R for n/ = 0; the relative viola- 
tions of rotational symmetry become smaller with bigger 
R and the one-loop violations are smaller than the tree- 
level ones. In the lower part of the figure the ratio of 
the two ratios is displayed, which amounts to replacing 



the continuum 1/R term that multiplies the logarithmic 
running of Vl,2 by a lattice "[1/R]" function. In doing so 
we isolate those lattice artefacts that appear specifically 
at order g 4 ; for instance, the lattice tadpole contribu- 
tions cancel from this combination. The sign is oppo- 
site to that of the tree-level differences, which explains 
the weaker relative violations at 0{g A ) in the comparison 
with 0(g 2 ) in the upper part of the figure. 




FIG. 11: The ratio of lattice and continuum order a and a 
pure gauge 5(7(3) potentials as a function of a 2 /r 2 for on-axis 
separations. 



C. The continuum limit 

The limit of large R = r/a has two interpretations: 
taken at a fixed lattice spacing a (sufficiently small for 
perturbation theory still to be reliable at a distance r) 
continuum pcrturbative predictions should be met. We 
have already discussed this scenario above and indeed 
demonstrated this agreement at large R in Figs. || and[ic|. 
On the other hand taking the limit of large R at fixed 
physical r corresponds to the continuum limit a — > 0. 
In Fig. [ll] we investigate the approach to the contin- 
uum limit by displaying ratios of lattice and continuum 
perturbation theory results Vi, i/V C i for i = 1,2 versus 
1/R 2 = a 2 /r 2 for pure SU(3) gauge theory and on-axis 
separations R > 2. For the Wilson gauge action we ex- 
pect lattice artefacts to be a polynomial in a 2 and indeed 
no linear term is found. The solid curves represent fits 
that are quadratic plus quartic in the lattice spacing a. 
For off- axis separations such as R || (1,1,0) we find a 
similar picture, but with different a 2 and a 4 coefficients. 

The same comparison is displayed for the fcrmionic 
contributions Vlj(R) — Vf(R) — V/(oo) and V c j(r) = 
limo-»o VL,f(r/a) alone for massless quarks in Fig. [l^. 
The fit curves are quadratic plus cubic for SW (csw — 1 
to this order in perturbation theory) and KS fermions 
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Fig. [L2| with Fig. [l3|, that corresponds to a quark mass 
ma = 0.1, reveals: "improvement" enables one to simu- 
late on coarser lattices at the same physical quark mass. 



0.05 0.1 0.15 0.2 0.25 0.3 

a 2 /r 2 

FIG. 12: Lattice artefacts on the fermionic contribution to 
the on-axis potential as a function of a 2 /r 2 for different ac- 
tions of massless quarks. Note that the leading order contri- 
bution is proportional to a/r for Wilson fermions. 
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FIG. 13: Lattice artefacts on the fermionic contribution to 
the on-axis potential as a function of a 2 /r 2 for ma — 0.1. 



while for Wilson fermions we indeed require the expected 
linear term and attempt a linear plus quadratic fit. The 
coefficient of the linear term turns out to be so small 
that the quadratic term already dominates for distances 
a/r > 0.1. We see that in spite of the more favourable 
functional dependence on a the perturbative lattice arte- 
facts of SW fermions are numerically bigger over the 
whole displayed R range than those of Wilson fermions. 
However, the KS action "outperforms" both Wilson and 
SW for this observable. As the quark mass is increased 
the violations of the continuum symmetry become more 
pronounced in the Wilson case while for the other two 
actions the change is only small, as a comparison of 
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FIG. 14: The same as Fig. ^| but for the fermionic contribu- 
tion alone. 
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FIG. 15: The same as Fig. but for ma = 0.1. 

From Fig. || it is obvious that for realistic values of nj 
the potential is dominated by the gluonic contributions 
and the additional violations of rotational symmetry due 
to sea quarks that we have just discussed are an interest- 
ing but numerically subleading effect. 

We shall now return to the limit where continuum per- 
turbation theory is met, i.e. we investigate the behaviour 
at fixed (small) a and large r = Ra. In the massless 
case this limit and the limit of large R = r/a at fixed 
r discussed above are equivalent, up to non-perturbative 
effects. However, as soon as a second scale ma > is 
introduced, the situation r > to" 1 becomes distinguish- 
able from r < to -1 . In Fig. [l4| we again display the 
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self energy subtracted fermionic contribution Vl,/(R.), 
this time multiplied by —R, in analogy to Fig. |[ This 
combination, multiplied by C^n/ isolates the shift that 
is induced onto an effective coupling Cf&r, due to sea 
quarks. The logarithmic slope of the three curves that 
are displayed in the figure is determined by the fermionic 
contribution to (3q and is therefore universal. The inclu- 
sion of sea quarks not only reduces /3 U5 relative to the 
pure gauge case, and hence slows down the running of 
an but it also decreases the absolute value of the 0(g 4 ) 
coefficient. In the case of massive quarks, however, this 
effect is (over) compensated by an increase in the coupling 
g 2 , relative to the pure gauge case, if we require the same 
physics at scales r ^> m~ x whe re the sea quarks decouple. 
This will be discussed in Sec. [V C 3 below. 



-0.001 - 



> 

QC 




FIG. 16: The same as Fig. [H] for Wilson fermions of various 



While for massless flavours the coupling runs loga- 
rithmically, for massive quarks the running gradually 
switches off around 3 R w 0.3 /(ma) w l/(-\/Uoma): 
physics at length scales Ra ^> mT 1 is insensitive to the 
presence of massive sea quarks, at least in perturbation 
theory. In Fig. |l^ we again compare the contributions to 
the effective Coulomb couplings resulting from the three 
different fermionic actions but for ma = 0.1. At this 
mass value Wilson and KS results accidentally happen 
to be very close to each other. In Fig. [l6] we investigate 
the mass dependence of the fermionic term somewhat 
more systematically by studying the example of Wilson 
quarks alone. In order to relate the lattice to a contin- 
uum scheme at large R, we not only have to subtract 
the (sea quark mass-dependent) static source self energy 
contribution to eliminate terms proportional to R from 
the figure but also the vertical offset, Ki(ma), at small 



R changes due to a mass dependent term that appears in 
the conversion factor b\ between the MS and the lattice 
scheme at fin ite a . We will discuss this effect in some 
detail in Sec. Ill D below. 

In the limit r 3> m^ 1 all curves indeed approach a con- 
stant value: the running of the effective coupling is not 
affected by the presence of massive sea quark flavours 
anymore. The resulting potential is the same as that in 
the pure gauge case, at least in perturbation theory, al- 
beit with a different overall normalization of the effective 
Coulomb coupling. We shall discuss this shift of the QCD 
coupling constant at a given scale and the possibility of 
matching to the quenched theory in Sec. Ill E below. 



While at large R the lattice data and the continuum 
curve nicely coincide with each other [once the offset 
K\ (ma) has been determined from the large R data and 
subtracted], at small R we observe large 0(am) lattice 
artefacts, in addition to mass independent 0(a) effects: 
even at the relatively small value ma — 0.03 the lattice 
points lie systematically below the large R expectation. 
Fig. [l5] confirms this qualitative pattern for the other two 
fermionic actions. 



D. 



"AATi" 



As we have seen above the choice of lattice action not 
only affects the quality of rotational symmetry at small 
R but also the overall normalization. The MS scheme is 
related to the lattice scheme via 

ttj^jfa -1 ) = (Xl + b\o? H , (40) 

with the conversion factor 

b 1 = -n/(2N) + k 1 N + K 1 (ma)n f . (41) 

The numerical constant k\ is known for a variety of glu- 
onic actions and K\(0) is known for Wilson, SW and 
KS q uark s a nd in dependent of the gluonic action, cf. 

Eqs. flAl6| ) - ( |A17| ). 

The coefficients of the /3-function in the MS'-scheme as 
well as in the lattice scheme in the continuum limit do not 
depend on the quark mass: both are "mass- independent" 
renormalization schemes. In lattice simulations it is of- 
ten worthwhile to analyse quantities prior to an extrap- 
olation to the continuum limit. One such example is 
determinations of the QCD running coupling from small 
Wilson loops jl], ||, ||, |j, H Hi that are ill-de fined in this 
limit. We will discuss this technique in Sec. IV B below. 



Another example is simulations within the framework of 
an effective field theory like NRQCD E2| that requires a 
finite momentum cut-off. 

At finite a and m > the lattice scheme becomes mass- 
dependent, as indicated by the argument of the function, 



Ki(ma) = Ki(0) + AK^ma), 



(42) 



The constant Co ~ 5.2 [H is introduced in Eq. (\3t) 



within Eq. (|4l]). In this case universality is lost and the 
coefficients of the perturbative /3-function acquire addi- 
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tional contributions that, like AK\, will in general de- 
pend on the dimensionful observable that is studied: 



1% 



MS 



dAK x (ma) 



(43) 



e?ln a 2 

By definition, Aifi(O) = 0. However, the behaviour of 
AKi is also constrained in the limit ma — > oo where the 
sea quarks decouple and therefore (3q' is = (11N)/(12tt) 
while $ = (Sjp 7 - n//(67r), i.e., 

AJfi m ^>°° - — \n(Dma) , (44) 

37T 

with a dimensionless constant D that we calculate in 



Sec. HIE below. 

We have seen above that in the limit r <C tvT 1 the be- 
haviour of the massless theory is emulated while the run- 
ning of the potential for r 3> m _1 is effectively quenched: 
the only scale that is relevant in perturbation theory in 
these limits is the distance r (or, equivalently, momen- 
tum q sw and therefore lattice artefacts disappear 
like (a/r) u [or (aq) v ] with some positive integer power v, 
whose value depends on the lattice action, as r — » oo. 
However, in the intermediate range of quark masses lat- 
tice corrections (am) 1 ' become relevant and the univer- 
sality of the /3-function is lost. 



The function K\(ma) of Eq. (A14) can in principle be 
read off from figures such as Figs. |15| and |l6| up to short 
distance lattice artefacts: Ki(ma) = — Cp/j An) Vl^D — 
y/C^ma + 0[(ma) 2 



,crj [cf. Eqs. (|37|), (|A5l|) and (|A76[) 1. 

For the (improved) SW action not only the lattice arte- 
facts are much more pronounced and qualitatively differ- 
ent from the Wilson and KS cases but also this overall off- 
set Ki is enhanced (although not its mass dependence). 



<3 
"CO 

< 




0.1 0.15 0.2 0.25 
ma 



FIG. 17: AK\ determined from the static position space 
potential, relative to Ki(0). 

We determine K\(ma) from two parameter fits, 



= Ao + ^ilnfi+^Ein 



[cy^maR 



TABLE X: Aifi as a function of ma for Wilson, SW and KS 
quarks. In the last row we display —K 1 (0) for comparison. 



ma 




AK\ ' 


AK[ ' 


0.01 


0.0005(13) 


0.0121(18) 


0.0038(20) 


0.03 


-0.0002(12) 


0.0316(15) 


0.0131(19) 


0.04 


0.0012(18) 


0.0394(18) 


0.0127(20) 


0.05 


-0.0006(27) 


0.0445(17) 


0.0079(26) 


0.10 


-0.0012(11) 


0.0632(17) 


0.0095(10) 


0.25 


0.0023(15) 


0.0837(14) 


0.0280(14) 


-Ki(Q) 


0.032983419 


0.0841444 


0.3957496 



n f AKi(ma) + 



(45) 



to the R > 9 on-axis data points, with one redundant 
parameter c whose role is to parametrize any residual 
lattice artefacts. The abov e fun ctio nal fo rm with the con- 
stant s Ap and A\ of Eqs. (A73) - (A75) is motivated by 



Eq. ( |A72| ) of Appendix]A~4 The special function Ein(x) 
deno tes t he normalized exponential integral defined in 
Eq. (|A52|) . 

All fits turned out to be stable within statistical er- 
rors under variations of the fit range as well as consistent 
with the off-axis data. The results on AK\ are displayed 
in Tab. [x|. In Fig. 17 these shifts are plotted relative to 
— ifi(Q), together with cubic splines. Since -Ki(O) < 0, at 
large ma the curves will diverge towards the negative di- 
rection. Note that the slope of AK\(ma) at ma = has 
also been determined in Ref. [0 for SW fermions, how- 
ever, by matching the vacuum polarization in momentum 
space, rather than the position space potential and AK\ 
is not universal. 

For the two 0(a) improved actions AK\ (ma) is small 
(a few %), relative to Ki(0), at least for ma < 1/4. This 
is very different in the case of Wilson fermions, where 
even at ma = 0.05 AK\ corrects K\(Q) by more than 
50 %. Using the static position space potential as the 
prescription that defines AA'i , the ma dependence is well 
parametrized by AK\(x) = ax + bx 2 with a = 1.20 ±0.02 
and b = — 5.7 ± 0.3 for x — ma < 0.1 in the Wilson case. 
For SW fermions the fit parameters read a = 0.44 ± 0.08 
and b = — 3.5±0.8 while the KS data are compatible with 
zero within our accuracy over the whole range x < 0.25 
and can be fitted by a straight line with a = — 0.04±0.04. 



E. The "/3-shift" 

We have seen above that at distances r ^> m~ x the 
running of the coupling is not affected by the presence 
of sea quarks anymore: at large distances, at least in 
perturbation theory, the effect of massive quarks can be 
integrated out into a shift of the coupling constant of the 
quenched theory. In contrast this is not possible for a 
theory with massless quarks, which completely decouples 
from the quenched case at all scales. 
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The theory with massive sea quarks and the quenched 
theory can be matched in the infra red by use of an 
intermediate mass-dependent scheme, such as the R- 
scheme, defined through the inter-quark potential in po- 
sition space, 



olr ( r 



(46) 



~^(R) (47) 

Inserting the perturbative expansion of Vint (R) hi terms 
of aL one obtains, 



„(°) , 



1 {R) a L 



(48) 
(49) 



where /i = (Ra)^ 1 , with known coefficients, e[ f \R). 
Note that these coefficients will depend on the function 
Ki(ma) discussed above. The requirement of a unique 
(static source self-energy subtracted) potential V ln t (R) at 
R ^S> (ma) -1 then results in, 



{l+ [e^ ) (oo)-e(° ) (oo)] a} a*£ f) + 



(50) 



One can now re-write the above equation in terms of 
the lattice parameter (3 = 2N/g 2 — N/(2kcil). For nj 
degenerate quark flavours with mass m and N = 3 the 
result reads, 

p{n f ) = pW +A p : ( 51 ) 

A/3 = ^r[\n(Dma) + 3TrAKi(ma)}, (52) 

Z7T 



F. Boosted perturbation theory and q* 

Lattice perturbation theory is well known for its bad 
convergence, partly due to large contributions from lat- 
tice tadpole diagrams. We recall that l^(R) incorporates 
such a contribution: (2A 2 -3)/(24A)T/i(R) [cf. Eq. ©]. 
Hence reordering the series in terms of a better behaved 
expansion parameter like av(q), the coupling defined by 
the static QCD potential in momentum space, 



V(q) 



-AnC F 



a v {q) 



(57) 



is desirable in many cases. In some respect this is similar 
to the situation in continuum perturbation theory and 
resembles an expansion in terms of a renormalized, rather 
than a "bare" , lattice coupling parameter. We will refer 
to such techniques as "boosted perturbation theory" . 




with the numerical values for the constant, 



3^(0)-- : (53) 
o 



D w = 0.448 ±0.002, (54) 
D sw = 0.0238 ± 0.0001, (55) 
D KS = 0.726 ±0.002, (56) 

for Wilson, SW and KS fermions, respectively. Since 
A(3 — > as ma —* oo, the ln(Dma) term has to be can- 
celled by 3nAKi (ma) at large ma in Eq. (|52|): the con- 



stant D is identical with that appearing within Eq. (44). 
We discuss the matching procedure in some more detail 
in Appendix [b|. 

Naively one would assume perturbation theory to be 
applicable as long as the relative a-shift is small. To 
leading non-trivial order in perturbation theory the f3- 
shift is linear in rif. At the next order, the situation will 
be complicated by additional te rms t hat are proportional 



FIG. 18: q*a values for the static potential. 

At 0(a 2 ) we can write, 

a v (q*) =a L + [a x +b x - 2/3 ln(q*a)\ a 2 , (58) 



with ai, b i and/3p as defined in Eqs. ( |A34[ ), ( |A3q ) , 
Eqs. ( |A1^ ) - ( |A16| ) and Eq. ( |A2| ), respectively. The 
"optimal" scale q* depends on the underlying process. 
It has been argued by Brodsky, Lepage and Macken- 
zie (BLM) J|, that the logarithmic average of the 
momenta exchanged at tree-level is a particularly good 
choice for the scale q* to be used within a one-loop calcu- 
lation. The scale optimization procedure can also be gen- 
eralized to higher order perturbative calculations [l5| . 
We illustrate the original recipe for the case of the po- 
sition space potential. Here the tree-level calculation, 
Eq. ©, yields, 



In [aq*(R)] 



to nj as can be seen from Eq. ( B20 ) of Appendix 



/^glngsin 2 ^) /q 2 
fd^qsm 2 (an*) /q2 



(59) 
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For R > 1 we find, 

*ru\ * i ^ 2 i * ln -^ i 
aq (R) = fx + — + h— + ■■■ 

with the numerical constants, 

/i w 1.672, / a w 1.00, / a « 0.41. 



(60) 



(61) 



The constants f 2 and have been obtained from a 
fit to R > 4 data while the r — > oo limit /i has 
been calculated directly. The fit curve and q* values 
are displayed in Fig. fTM Note that violations of rota- 
tional invariance in <7*(R) are remarkably small. Anal- 
ogously we can obtain q* values for the inter-quark 
force. We display results for the plaquette □, the chair, 
the "parallelogram", the "rectangle" VF(1,2), the static 
self energy <5m s t a t = Vs/2, the potential and the force 
F RlR2 = a[V(R^)- V{R 1 )]/(R 2 - Ri) at selected dis- 
tances in Tab. KH In addition, the q* values for the 
"tadpole improved" quantities m sta t + In □/(4a) as well 
as V(l) + lnD/(2a) and V{2) + lnD/(2a) are included. 

One can further convert between the V and the MS 
scheme: 

"ms(m) = «v(9*) - [oi + 2/3 ln(/x/«*)] a 2 + ■ ■ ■ . (62) 
It has been argued |43) that the scale 



ii = e- 5 / 6 q* 



(63) 



at which the n / dependence of a\ above is cancelled (in 
the massless case) was the optimal choice of scale for this 
conversion. Note that the ratio /i/q* above is indepen- 
dent of the number of colours N or flavours rif. The /i 
values are also included into the table. 

The average plaquette turns out to be the most ul- 
tra violet quantity, with q* ss 3.402/a, followed by the 
potential V(R = 1). Due to the self energy contribu- 
tion Vs = V(oo) = 2<5m s tat, q* is bounded from below 
by 1.672/a for the potential at all distances. However, 
in the case of the force q* F — > as R — > oo. As one 
would expect q* also approaches zero at large distances 
for the potential if Vs is subtracted. We find for instance 
q*a k, 0.59 for V(l) - V s as opposed to q*a w 2.86 for 
V(l) alone. 



G. The static source self energy 



We define, 



V s = lim V(R) 

R—>oc 



2Sr 



(64) 



where <5m s t a t is the lattice pole mass of a fundamental 
static colour source, e.g. a quark in the m — > oo limit. 
The above relation only holds in perturbation theory 
where the interaction energy vanishes as r — > oo. 

The static (or residual) mass <5m s t a t has been calcu- 
lated to 0(a 3 ) for pure gauge theory [Ell Q as well as to 



TABLE XI: 


BLM scales q* for some quantities. 


quantity 


q* a 


fj,a = e 5 / 6 g*a 


□ 


3.402 


1.478 


chair 


3.300 


1.434 


"parallelogram" 


3.128 


1.360 


W(l,2) 


3.066 


1.332 


<5m s tat 


1.672 


0.727 


V(l) 


2.860 


1.243 


V(2) 


2.317 


1.007 


Fl2 


1.025 


0.445 


F 23 


0.904 


0.393 


Srristat + lnD/(4o 


0.835 


0.363 


V(l) + lnD/(2o) 


1.701 


0.739 


V(2) +lnD/(2o) 


1.316 


0.572 


V(V2) 


2.558 


1.112 


v(Vs) 


2.417 


1.050 


V(3) 


2.173 


0.944 


V(4) 


2.062 


0.896 


V(5) 


2.007 


0.872 


V(6) 


1.959 


0.851 


V(7) 


1.930 


0.839 


V(8) 


1.902 


0.826 


1/(16) 


1.805 


0.784 


V(24) 


1.768 


0.768 


V(28) 


1.757 


0.764 


V(32) 


1.748 


0.760 



0(a 2 ) for massless Wilson and SW quarks ||. This has 
enabled a number of authors to obtain the b quark mass 
mt,(rrTb) from lattice simulations of static-light mesons. 
Since the sources are static the 0(a) result can be ob- 
tained from a tree-level calculation and we count this 
as leading order (LO), discarding the trivial 0(1) value 
Srristat = 0. The one-loop 0(a 2 ) result is then next 
to leading order (NLO) and the two-loop 0(a 3 ) value 
NNLO. Note that the counting conventions employed in 
some of the literature differ from the one defined above, in 
which 0(a n ) corresponds to an N m LO calculation with 
m = n — 1 for which diagrams involving m loops have to 
be computed. For instance in Ref. S the value m = n 
is used, creating the (wrong) impression that a two-loop 
calculation is required to obtain the 0(a 2 ) result. 

We have computed the static quark mass shift as a 
function of the mass of three species of sea quarks: KS, 
Wilson and SW. The results, most of which can also be 



read off from the last seven lines of Tab. |IX[ are displayed 
in Tab. XII. The labelling conventions are identical to 
those of Eqs. @ - ©: 



aSm st at = SM ig 2 + SM 2 g 4 , (65) 
SM 2 = 8M pg + C F n f 5M f , (66) 
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TABLE XII: Fermionic contributions to the static self energy 
of a fundamental source. The labelling conventions are de- 
fined in Eqs. @ - (§§). 



ma 


oivif /ru 


OJVLf /tU 


OJVLf /tU 


olvlf /tU 





-1.1667 (2) 


-1.3404 (2) 


0.1633 (5) 


-1.8587 (6) 


0.01 


-1.500 (3) 


-1.3061 (2) 


0.1489 (1) 


-1.8540(11) 


0.03 


-1.142 (3) 


-1.2402 (2) 


0.1216 (2) 


-1.8443 (2) 


0.04 


-1.128 (2) 


-1.2090 (1) 


0.1083 (6) 


-1.8397 (5) 


0.05 


-1.120 (3) 


-1.1791 (1) 


0.0967 (5) 


-1.8364 (1) 


0.10 


-1.0801 (1) 


-1.0447 (1) 


0.0478 (1) 


-1.8134 (1) 


0.25 


-0.922 (2) 


-0.7541 (1) 


-0.0268(1) 


-1.7364 (1) 


1 


-0.2962 (1) 


-0.2336 (2) 






2 


-0.06131(3) 


-0.08264(6) 






4 




-0.02030(2) 







CO 

a? 




0.25 





r 




-0.2 




-0.4 - 


CO 

O 


-0.6 - 


to 


-0.8 - 




-1 - 




-1.2 - 




-1.4 - 




FIG. 19: The fermionic contribution to the static self energy 
Sm stat for Wilson and KS quarks as a function of the quark 
mass. 



where 



8M X 



C, 



Moo), 



N 



SM pg = i- [Vh(oo) + V NA (oc)} 



2 ^ 2 

2N 2 -3 



247V 



V T (oo) 



(67) 



(68) 



and SMf = SM\ K ' for KS quarks and 



SM f = 5Mf ] + cswSMy + c 2 sw SM { /> (69) 
for Wilson- SW fermions. We find the numerical values, 



SMx w 0.126365504 C F 
5M pq w 0.05292 C F , 



(70) 
(71) 



FIG. 20: The fermionic contribution to the static self energy 
for Wilson, SW and KS quarks, normalized to the respective 
zero quark mass results. 



where the latter numbe r ap plies to SU (3) gauge theory 
(cf. the last line of Tab. [Vl]). Note that this result differs 
in the least significant digit from the value obtained in 
Ref. that can be translated into 5M pg w 0.05297C F . 
The mass dependence of SMf [which in Eq. (^6|) is ac- 
companied by a factor C F nf] is visualized in Fig. [l9] for 
Wilson and KS quarks. SM f approaches zero like l/(ma) 
as ma — > oo. 

Finally, in Fig. |2(] we compare the results on SMf (ma) 
for quark masses ma < 0.25, normalized to SMf (0), from 
all three quark actions with each other. While in absolute 
terms SMf for SW quarks turns out to be much bigger 
than for the other two actions (cf. Tab. XII) we find the 
relative variation with the quark mass to be much weaker 
for both 0(a) improved actions than for Wilson fermions. 
This is consistent with our observations for small Wilson 
loops and A_?Ci above. For ma < 0.1 5Mf(ma) is well 
parametrized by a quadratic function, 



SMf (ma) = SMf(0) + bma + c(ma) 



(72) 



with SMf(0) = -0.0011667(3), -0.0013404(2) and 
-0.0030351(11), b = 0.00101(8), 0.00350(1) and 
0.00245(5) and c = -0.0014(8), -0.00543(5) and 
—0.0020(4) for KS, Wilson and SW quarks, respectively. 
Starting from the expansion of the plaquette, 



□ = 1 — ClQfL 



c 2 a 



(73) 



we can rearrange the above expansion of the self energy 
in the following way, 



aSm s 



lnD 
lnD 



m Y a L 



m\a\ 



(74) 
(75) 



17 



TABLE XIII: Fermionic contributions to the static self energy 
as defined in Eqs. ( ]74| ) - (Jts|) . The first three rows refer to 
the expansion of the tadpole improved observable in terms of 
cxl, the last rows to an expansion in terms of ajjg-(0.363/a). 





KS 


Wilson 


sw 


a f 


-0.1853 


-0.2090 


-0.5051 


bf 


0.21 


0.61 


0.39 


c f 


-0.38 


-1.02 


-0.33 


0/ 


-0.03499 


-0.05860 


-0.3548 


bf 


0.26 


-0.67 


-0.081 


Cf 


-0.38 


5.08 


3.41 



where /! = q* 



a—5/6 



0.363/a, m\ = (4tt)<5Mi - c x /4, 



m\ = {AttY8M 2 - c 2 /4 - cf 



We will refer to prac- 



tices such as adding the In □ term non-perturbatively 
and subtracting it perturbatively as "tadpole improve- 
ment" of the observable while in the second step, Eq. (|75| ) 
we "boost" the perturbation theory. Here we choose to 
convert everything into the MS scheme. We prefer this 
to plaquette based schemes or the V scheme as now all 
mass dependence is made explicit in the coefficient func- 
tion m^ma). A conversion of the result into another 
scheme of choice including mass dependent schemes is 
easily possible with the help of Appendix A]. We will 
obtain the MS couplin g from the average plaquette and 
other quantities in Sec. IV B below. 



For SU(3) gauge theory with rif quark flavours of 
masses ma < 0.1 we find, 

m\ = 1.0700768, (76) 
ra\ = 7.6104 + n } [a } + b } ma + c f (ma) 2 ] , (77) 
ml = -0.5839 + n/ [a/ + bfma + c f (ma) 2 ] , (78) 



of Tab. pffll 



Tadpole 
11.14 



with the constants aj, bf, Cf, 
improvement reduces the pure gauge value 5m pg 
to TOp g « 7.61 and boosting reduces this further down to 
rfipg « —0.584. The fermionic coefficients 5m f i— > a/ i— > 
5/ also undergo a reduction in each step. The NNLO 



coefficient Sm^ « 86 



67.5 



2.85 is only 



known for n 



f 



0. 



When applying the above result to extract say mj from 
lattice simulations in the static limit one can readily 
ignore the mass dependent terms for Wilson fermions, 
which are 0(a) effects. However, it does not harm to 
include them either. In the cases of massive KS and SW 
sea quarks which are 0(a) improved at least the 6/s have 
to be included. 



IV. COMPARISON WITH 
NON-PERTURBATIVE DATA 

We shall compare and apply our perturbative calcula- 
tions to non-perturbative results obtained in lattice sim- 
ulations. For this purpose we will use several data sets 




that have been obtained by different collaborations. All 
quenched reference data are from simulations of one of 
the co-authors of this article and collaborators (Ref. E3] 
and references therein) while rif = 2 data have been 
provided by the SESAM/TvT Collaboration 
(Wilson fermions), the UKQCD Collaboration 
(SW fermions) and the MILC Collaboration 
fermions 4 ). 

We will apply our perturbative results to the "/3-shift" 
encountered when including sea quarks — relative to the 
quenched approximation. Subsequently we determine 
the QCD running coupling "constant" for nf = and 
rif = 2 from lattice data of the correlation length ^] 
Ro = r o I a , the average plaquette and the short distance 
static potential. In this context we will also make use of 
the CP-PACS @ ensemble obtained with n f = 2 SW 
fermions and the Iwasaki gluonic action [ p6[ . Finally we 
shall also compare perturbative and non-perturbative lat- 
tice potentials with the aim of parametrizing lattice arte- 
facts and to resolve the differences in the running of the 
Coulomb couplings between quenched and un-quenched 
data sets. 



A. The non-perturbative "/3-shift" 

We study the situation of n/ mass-degenerate flavours 
of sea quarks. Two parameters can be varied: the cou- 
pling j3 = 6/g 2 and the lattice quark mass ma which in 
the case of Wilson-SW fermions is related to the parame- 
ter k. Each k-(3 (or ma-(3) combination can be translated 
into a pion mass m T ro and a lattice spacing a/ro where 
ro is a correlation length defined below. As — * oo we 
reach the continuum limit, a/ro — * 0. Sending the quark 
mass mro to zero corresponds to a vanishing pion mass, 
?7i 7*0 — > (at finite (3: only up to violations of chiral sym- 
metry). However, in general the two limits will "mix": 
varying [3 at fixed ma will not only affect a/ro but also 
to some extent the ratio m^r® while a variation of ma, 
keeping the coupling fixed, will result in a change of a/ro 
as well. 

In view of the computational cost of lattice simulations 
incorporating sea quarks, predicting by what amount the 
coupling has to be shifted in order to compensate for the 
change in a/ro that is induced by varying the quark mass 



4 While KS quarks are only defined for multiples of four mass- 
degenerate fermion flavours, many authors have attempted to 
emulate nj = 2 (or even rif = 3) by using the (positive) square 
root of the rif = 4 fermionic determinant in their simulation. 
The MILC Collaboration adopt this strategy. However, it is not 
clear that the resulting lattice action corresponds to a local field 
theory like QCD. In our analysis we shall set this problem aside. 
We remark that the perturbation theory generated by the use of 
this action indeed corresponds to replacing nj = 4 by tij^ = 2 
to all orders, at least as long as no external fermion lines are 
encountered. The comparison with simulation results, however, 
might be meaningless should no universal continuum limit exist. 
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ma is certainly desirable. We will explore this possibility 
by comparing the expectation Eq. (q2|) on the perturba- 
tive /3-shift against results from lattice simulations with 
rif = 2 sea quarks of masses 0.0098 < ma < 0.3. 

At large distances the potentials will be dominated by 
non-perturbative effects: the quenched potential will lin- 
early rise ad infinitum. However, as soon as sea quarks 
are introduced the Z 3 symmetry of the action is broken 
and at some distance |39| r c w (2.10 + 1.56mro)fo string 
breaking will set in; obviously the linear quenched be- 
haviour at large distances can never be emulated by a 
theory with a completely flat potential at large r > r c . 
Sea quarks will also affect the running of the coupling 
at short distances r < 0.3 m _1 . Hence, the best we can 
hope for is that sea quarks decouple within a window of 
distances m _1 <C r < r c . 



TABLE XIV: Comparison between non-perturbative /3-shifts 
A/3 ro and the respective perturbative predictions A/?n for 
two flavours of Wilson fermions [B9 , ra . 



/3 


ft 


ma 


ro/a 


-A/3 ro 


-Aft., 


5.5 


0.158 


0.0597 


4.03(3) 


0.352(5) 


0.318(3) 


5.5 


0.159 


0.0398 


4.39(3) 


0.395(5) 


0.371(4) 


5.5 


0.1596 


0.0280 


4.68(3) 


0.427(4) 


0.416(4) 


5.5 


0.160 


0.0202 


4.89(3) 


0.450(4) 


0.456(4) 


5.6 


0.156 


0.0504 


5.10(3) 


0.373(5) 


0.342(3) 


5.6 


0.1565 


0.0402 


5.28(5) 


0.399(5) 


0.369(4) 


5.6 


0.157 


0.0300 


5.46(5) 


0.410(5) 


0.406(4) 


5.6 


0.1575 


0.0199 


5.89(3) 


0.454(5) 


0.457(4) 


5.6 


0.158 


0.0098 


6.23(6) 


0.488(7) 


0.538(5) 



TABLE XV: Comparison between non- 
A/3,- and the respective perturbative 
two flavours of SW fermions Effl. 



-perturbative /3-shifts 
predictions A/3i_i for 



/3 


K 


ma 


ro/a 


-A/3 ro 


-A/3!., 


5.2 


0.135 


0.0459(2) 


4.75(4) 


0.736(5) 


0.679(3) 


5.2 


0.1355 


0.0236(2) 


5.04(4) 


0.766(5) 


0.750(3) 


5.2 


0.13565 


0.0179(9) 


5.21(5) 


0.784(6) 


0.780(6) 


5.25 


0.1352 


0.0427(2) 


5.14(5) 


0.723(6) 


0.686(3) 


5.26 


0.1345 


0.0720(2) 


4.71(5) 


0.671(5) 


0.632(4) 


5.29 


0.134 


0.0927(3) 


4.81(5) 


0.641(5) 


0.609(3) 


5.29 


0.1350 


0.0535(2) 


5.26(7) 


0.699(7) 


0.664(3) 


5.29 


0.1355 


0.0350(1) 


5.62(9) 


0.736(9) 


0.708(3) 



We shall define a 

Ap ro (n f ,r /a,ma) = ^ n f\r a /a,ma) - /3 {0) (r /a), 

(79) 

by non-perturbatively matching un-quenched and 
quenched (3 values that correspond to the same scale fill 



TABLE XVI: Comparison between non-perturbative /3-shifts 
A/3 ro and the respective perturbative predictions A/3i_i for 
"two" flavours of KS fermions @. The /3 = 5.6, ma = 0.025 
data point is from Ref. |52j. 



/3 


ma 


ro/a 


-A/3 ro 


— Ajfli-i 


5.3 


0.3 


1.65(3) 


0.079(25) 


0.149(3) 


5.35 


0.3 


1.79(1) 


0.084(20) 


0.149(3) 


5.415 


0.3 


1.97(5) 


0.076(9) 


0.149(3) 


5.3 


0.2 


1.75(5) 


0.119(25) 


0.196(2) 


5.35 


0.2 


1.87(1) 


0.111(13) 


0.196(2) 


5.415 


0.2 


2.15(1) 


0.125(10) 


0.196(2) 


5.3 


0.15 


1.78(1) 


0.130(20) 


0.226(3) 


5.35 


0.15 


1.94(1) 


0.132(12) 


0.226(3) 


5.415 


0.15 


2.27(1) 


0.154(11) 


0.226(3) 


5.3 


0.1 


1.86(2) 


0.157(15) 


0.267(2) 


5.35 


0.1 


2.04(1) 


0.161(10) 


0.267(2) 


5.415 


0.1 


2.45(1) 


0.194(11) 


0.267(2) 


5.5 


0.1 


3.10(4) 


0.225(7) 


0.267(2) 


5.3 


0.075 


1.91(2) 


0.172(15) 


0.295(2) 


5.35 


0.075 


2.15(1) 


0.190(10) 


0.295(2) 


5.3 


0.05 


1.96(2) 


0.188(13) 


0.336(2) 


5.35 


0.05 


2.35(3) 


0.237(13) 


0.336(2) 


5.415 


0.05 


2.72(1) 


0.246(9) 


0.336(2) 


5.5 


0.05 


3.42(3) 


0.273(6) 


0.336(2) 


5.3 


0.025 


2.03(1) 


0.209(11) 


0.406(1) 


5.35 


0.025 


2.36(1) 


0.240(12) 


0.406(1) 


5.415 


0.025 


2.95(1) 


0.286(8) 


0.406(1) 


5.5 


0.025 


3.80(2) 


0.324(5) 


0.406(1) 


5.6 


0.025 


4.8(1) 


0.341(11) 


0.406(1) 


5.415 


0.0125 


3.12(2) 


0.306(7) 


0.476(1) 


5.5 


0.0125 


3.98(1) 


0.346(4) 


0.476(1) 


5.6 


0.08 


4.08(1) 


0.259(4) 


0.289(4) 


5.6 


0.04 


4.54(2) 


0.312(4) 


0.358(1) 


5.6 


0.02 


4.84(1) 


0.345(3) 


0.429(1) 


5.6 


0.01 


4.99(1) 


0.361(3) 


0.499(1) 



0.5 fm, implicitly defined through the relation, 
odV(r) 



ro 



= 1.65. 



(80) 



If non-perturbative effects cancel from Eq. ( |79[ ) and lead- 
ing order perturbation theory applies then A(3 ro will only 
depend on ma (and nf) but not on the lattice spacing 
a/rQ. Perturbation theory will obviously break down for 
ma <C 1 in which case the coefficients of the perturbative 
expansion explode [cf. Eq. fl52|)]. This is a reflection of 
the fact that the running of the coupling in the theory 
with massless sea quarks differs from the nf — case at 
all scales. 

On the other hand, unless ma <C l/D the relative 
importance of the non- universal term AKi(ma) will in- 
crease and in general the matching of the running of 
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FIG. 21: The non-perturbative /3-shift for two flavours of 
Wilson fermions, in comparison with the perturbative predic- 
tion (curve). 
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FIG. 22: The non-perturbative /3-shift for two flavours of 
SW fermions. 



the inter-quark force at a scale ro is not equivalent to 
the matching of the running of the potential at larger 
distances anymore. In principle one can work out the 
matching condition for ro in perturbation theory. It 
turns out, however, that at distances ro w 0.5 fm the 
non-perturbative contribution to the lattice potential is 
already substantial. For instance a perturbative nt — 
calculation yields ro/a — 5.33 at (3 ps 1.49 while non- 
perturbatively this result is obtained at (3 ps 6.0 and (3 
depends only logarithmically on a! 

Given the above situation, it is difficult to iden- 
tify a sensible way of combining perturbative and non- 
perturbative results at large distances. Consequently, we 
refrain from attempting to do this but separately employ 
purely perturbative and purely non-perturbative match- 



FIG. 23: The non-perturbative /3-shift for "two" flavours of 
KS fermions. 



ing strategies: on the lattice simulation side of things we 
match r /a. In perturbation theory we match a quan- 
tity that is suitable for perturbative treatment, namely 
the inter-quark potential at large distances, rather than 
a quantity that is inspired by (non-perturbative) phe- 
nomenology, like ro- We then hope that some, ideally 
most, non-perturbative effects cancel each other on the 
right hand side of Eq. (jfjj) and leave the /3-shift within 
the available window of ro/a and ma untouched. Clearly, 
the perturbative matching will break down if either ma 
or f3 are too small. In addition the strategy that we adopt 
requires ma < 1 / D to limit violations of universality as 
well as mro > 0.3 to ensure that the sea quarks do not 
affect the running of the coupling at distances around ro- 
The quenched (3 value that corresponds to a given ro /a 
is determined by use of the interpolation, 



ro/a = exp [do + d\x + d 2 x 2 + d%x 3 ) 



5.9, 
(81) 

with d ps 1.489, d x w 1.982, d 2 ~ -0.630, d 3 pa -1.522, 
obtained from a fit to lattice results within the win- 
dow jlo| 5.5 < /? < 6.2. In principle we can also use 
the more recent precision data of Refs. p0|, 51 . How- 



ever, the difference is insignificant in view of the present 
level of accuracy of un-quenched data and the match- 
ing of some of the MILC data requires an interpolation 
that extends to lattice spacings coarser than those inves- 
tigated in the latter two references. In the un-quenched 
simulations, different groups used different procedures to 
extract ro/a which partly explains why different collab- 
orations can obtain very different error estimates with 
similar computational efforts. Aiming only at a quali- 
tative comparison we shall not attempt to reanalyse all 
data in one and the same way but prefer to cite the pub- 
lished values and errors. 

In the case of KS fermions ^8|, |32| there is no additive 
mass renormalization. For Wilson fermions lB9l, Em we 
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define the lattice quark mass via, 
1 f 1 

ma= 2{- K -; 



(82) 



where n c has been obtained from a chiral extrapolation 
of the non-perturbative m\ at fixed 5 /3 6 . In the case of 
the SW data §, £7) we do not know n c but we can use 
the published values of the PC AC quark mass [^|, which 
is equivalent to any other definition to leading order per- 
turbation theory. 

We compile results from the three collaborations in 



Tabs. XIV 



KVI 



We have sorted Tab. |XVI| (with 
the exception of (3 = 5.6) with respect to the quark 
mass and then (3 while the other tables are sorted 
the other way around. Unfortunately, no PCAC mass 
value is available for the most critical UKQCD data set 
[(13, k) = (5.2,0.13565)], however, we arrived at the es- 
timate ma = 0.0179(9) by extrapolating ma as a linear 
function of (m n ro) 2 to the value |03| ( m 7r ro) 2 = 0.132(2). 
The perturbative prediction Eq. (p2[) is included in the 
last column of the tables where the error is due to the 
uncertainty in AA'i(mo). 

In Figs. [2l| - ^ we compare the non-perturbative re- 
sults with the expectations. In all cases ma is much 
smaller than 1/D as desired. However, mro, varying 
from 0.06 to 0.24 in the Wilson case, from 0.09 to 0.46 
in the SW case and from 0.04 to 0.6 in the KS case, is 
not exactly big relative to r$ = 0.3 to -1 (the distance 
below which the perturbative running of the coupling 
changes and sea quark effects cannot completely be com- 
pensated for by a scale redefinition alone). Nonethe- 
less, the first two figures reveal an excellent agreement 
with the expectation within 10 % for Wilson-type quarks, 
even for such light quark masses. Moreover, no signifi- 
cant /3-dependence is observed, again in agreement with 
Eq. (p2|), indicating higher order corrections to be small. 
This is very different for the KS fermions depicted in 



5 In principle other extrapolations are possible like keeping r^/a 
fixed Wm- However, the whole matching idea is based on a semi- 
quenenca philosophy: what value of the rif = coupling will 
produce the same infra red physics, e.g. ro/a? Since at k c the 
low energy physics will be very different anyway, the "natural" 
choice in this case seems to be keeping the coupling constant 
fixed. Having said this, to the order of perturbation theory at 
which we work both approaches are equivalent anyway. 

6 To all orders in perturbation theory ma is defined via Eq. (p^). 
To the order that we work at, re c , P crt = 1/8, which means that 
the k values employed in the lattice simulation all correspond to 
negative masses, ma, when directly plugged into the perturbative 
expansion. Obviously this is not a sensible choice, and hence we 
formulate our perturbation theory in terms of the mass ma rather 
than k. Subsequently, we determine the n that corresponds to a 
given ma value non-perturbatively. In the case of SW fermions 
we proceed in an analogous way: on the perturbation theory side 
we use ma and csw = 1, which is the value consistent with the 
order at which we work. However, in the simulation we use the k 
and csw values respectively which result in the same ma value 
and that eliminate O(a) lattice artefacts non-perturbatively. 



Fig. |23|: here the agreement with the prediction only im- 
proves as the lattice spacing is reduced. 
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FIG. 24: /3-shift: the non-perturbative result, normalized to 
the prediction. 

In Fig. ^ we finally compare the three formulations 
within the window 0.1 > ma > 0.009 for 4 < r^ja < 6.5 
data. The Wilson and SW results seem to fall onto al- 
most universal curves that differ by less than 10 % from 
the prediction while the KS results deviate by much more 
and some /3-dependence is evident, Whether this is due 
to a slower convergence of the perturbative series, due to 
the inexact updating algorithm employed or due to nf 
not being a multiple of four is an open question. Around 



0.05 rn 



20 MeV the reliability of the match- 



ing for Wilson-SW quarks finally appears to break down 
(left- most data point in Fig. p4| ) : the behaviour becomes 
"truly un-quenched" . 

We remind the reader that /3-shifts are in general 
not independent of the quantity that is used to match 
quenched and un-quenched theories. This ambiguity ex- 
ists non-perturbatively as well as in perturbation theory. 
The qualitative agreement between prediction and sim- 
ulation for Wilson- type quarks indicates that, at least 
at present masses, physics at hadronic scales is not yet 
strongly affected by quark loops, which is consistent with 
the phenomenological success of the quenched approxi- 
mation. 

In simulations with non-perturbatively improved SW 
quarks the lines of constant a/ro and constant m^o & r e 
significantl y t ilted with respect to both axes of the f3- 
ma plane [p4| . This observation is consistent with the 
small value, D sw » 0.0238, but causes practical prob- 
lems as going to lighter quark masses at sensible a/ro 
values requires simulations at small (3s for which the 
non-perturbative determination of the improvement co- 
efficient csw causes problems plfl . In the worst case it 
might even be conceivable that the slope of the variation 
of a/ro with k — > k c at fixed (3 eventually diverges and 
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that the continuum limit [3 — > oo does not at all exist for 
light quarks |^6|, |5?j . The variation of r with the quark 
mass, however, is reduced when actions with only mod- 
erately negative values of Ki(0) [and therefore D values 
that are of order one] like the Wilson or KS action are 
employed. 

Our result might be taken as an indication that the 
perturbative matching of rif — 1 to n/ QCD couplings 
in the MS scheme via the intermediate mass-dependent 
F-scheme [[l4| (or by calculating and matching a physi- 
cal amplitude |58| ) is likely to be quite reliable, even at 
masses as light as or lighter than that of the charm quark. 
In this case the matching condition reads, 



1 



1 



i 

3^ 



hi 



'C - 



which may be rewritten as, 
where 



MS lw MS 



(M), 



(83) 
(84) 



number of quark flavours and masses is simulated the re- 
sulting a should become independent of the choice of the 
experimental input quantity. In this way the strong cou- 
pling constant can be determined from low energy hadron 
phenomenology. 

There arc of course higher order perturbative as well 
as non-perturbative corrections to Eq. ( |86| ) which, how- 
ever, will vanish as a — > 0. In practice these corrections 
are big in the range of lattice spacings that will ever 
be realistically accessible || |5j|. One way to improve 
the convergence is to convert between the MS and the 
lattice schemes at an optimized BLM scale ^|, [h| H 
fx = e~ 5 / 6 g*, rather than at a" 1 . We refer to this reorga- 
nization of the perturbative series as boosted perturba- 
tion theory. Another possibility (that can be combined 
with the BLM scheme) is to "measure" the coupling on 
the lattice. This can either be done from a quantity that 
depends on the lattice spacing like the average plaque- 
tte [Q, ||, ||, 0, H or at a scale fj, <C a -1 from quantities 
that have a well defined continuum limit |^9|, [H], ^2| . 
Here we shall follow the former strategy. 



p, = VCocxp [ -- ] m= 0.990(3)m, 



(85) 



to this order in perturbation theory 7 . Clearly, the appli- 
cability of the matching method in this case is ultimately 
limited by the reliability of the perturbative running of 
the MS coupling at small scales. 



B. Determining otjjg 

Lattice simulations yield correlation lengths and 
masses that are functions of a set of input parameters, 
namely the inverse lattice coupling (3 and lattice quark 
masses ma (or in the case of Wilson-type fermions re- 
values that can be related to the quark masses). One can 
then for instance obtain the strong coupling constant in 
the MS scheme, 



a- 



MS 



(a 



OLL 



(86) 



once a physical scale has been assigned to the lattice 
spacing a. One such possibility would be to "measure" 
Ra = ro/ a on the lattice and to equate rg ss 0.5 fm. 
Other input scales with a more direct connection to ex- 
periment are possible, for instance the proton mass tun 
or the pion decay constant f^. If QCD with the right 



Eqs. Q - Q can also be cast into a< n /'(m) = a( n / -1 >(m) + 
co 2 where tnc numerical value c £S 0.0033 is of the same mag- 
nitude as p3 c = 7/(727r 2 ) £3 0.0099. Note that the difference 
between pole and MS masses is irrelevant to this order in per- 
turbation theory. 



TABLE XVII: Pure gauge SU(3) results on ro and the average 
plaquettc □. 



p 




ro/a 




□ 


5.5 




2.01(3) 




0.49680(2) 


5.6 




2.44(6) 




0.52451(3) 


5.7 




2.86(5) 




0.54919(3) 


5.8 




3.64(5) 




0.56765(2) 


5.9 




4.60(9) 




0.58184(2) 


6.0 




5.33(3) 




0.59368(1) 


6.2 




7.29(4) 




0.61363(1) 


6.3 




8.39(7) 




0.62243(1) 


6.4 




9.89(16) 




0.63064(1) 


6.6 




12.73(14) 




0.64567(1) 


TABLE XVIII: Pure 


gauge 5(7(3) results on 


the short dis- 


tance potential and force J ? fl 1 H 2 


= a[V(R 2 ) - 


V(Ri)]. 


P 


aV(l) 


aV(2) 


F 12 


F 23 


5.5 


0.636 (2) 


1.155(21) 


0.519(22) 


0.29 (17) 


5.6 


0.566 (1) 


0.963 (4) 


0.398 (5) 


0.325(16) 


5.7 


0.506 (2) 


0.804 (5) 


0.297 (5) 


0.233(15) 


5.8 


0.464 (1) 


0.706 (3) 


0.242 (3) 


0.160 (6) 


5.9 


0.4330(9) 


0.645 (2) 


0.212 (2) 


0.127 (4) 


6.0 


0.4111(2) 


0.5974(2) 


0.1863(2) 


0.1032(4) 


6.2 


0.3778(1) 


0.5337(1) 


0.1559(1) 


0.0783(2) 


6.4 


0.3514(2) 


0.4889(4) 


0.1375(4) 


0.0648(5) 


6.6 


0.3293(2) 


0.4538(3) 


0.1244(3) 


0.0546(5) 
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TABLE XIX: Average plaquette, short range potential and 
force with two flavours of Wilson quarks. The corresponding 
ro/a and ma values can be found in Tab. XIV. 



TABLE XXI: The same as Tab. XIX but with "two" flavours 



of K S qua rks. The corresponding ro/a values can be found in 



Tab. 



p 


K 


□ 


aV(i) 


aV(2) 


Fl2 


5.5 


0.158 


0.55547(4) 


0.4756(3) 


0.7156(10) 


0.2400(10) 


5.5 


0.159 


0.55815(3) 


0.4685(3) 


0.6965 (9) 


0.2280 (9) 


5.5 


0.1596 


0.55967(3) 


0.4644(2) 


0.6853 (7) 


0.2209 (7) 


5.5 


0.160 


0.56077(5) 


0.4616(3) 


0.6761 (8) 


0.2145 (8) 


5.6 


0.156 


0.56989(1) 


0.4460(2) 


0.6511 (6) 


0.2051 (6) 


5.6 


0.1565 


0.57073(1) 


0.4440(3) 


0.6450 (9) 


0.2010 (9) 


5.6 


0.157 


0.57160(1) 


0.4422(3) 


0.6387(11) 


0.1965(10) 


5.6 


0.1575 


0.57257(1) 


0.4394(2) 


0.6336 (6) 


0.1942 (6) 


5.6 


0.158 


0.57337(1) 


0.4373(3) 


0.6303(13) 


0.1930(12) 



1. The method 

We collect the average plaquette, the potential at R = 
1 and at R — 2 obtained in quenched simulations as well 



as from simulations with sea quarks in Tabs. XVII 
In addition we include the "force" , 

F 12 = a[V{2)-V(l)], 
F 23 = a[V(3)-V{2)]. 



XXI 



(87) 
(88) 



We will restrict our discussion to a one-loop determi- 
nation of the running coupling. We shall also calculate 
two-loop corrections for the pure gauge case. Two-loop 
results are also known for the case of the plaquette with 
massive Wilson quarks 8 p9|. 



TABLE XX: The same as Tab. |XIX but with two flavours of 
SW quarks. The corresponding rp /a and ma values can be 
found in Tab. |XV. 



p 


K 


□ 


aV(l) 


aV{2) 


F 12 


5.2 


0.135 


0.53368(1) 


0.4823(2) 


0.6970 (8) 


0.2147 (8) 


5.2 


0.1355 


0.53629(1) 


0.4762(2) 


0.6832 (8) 


0.2070 (8) 


5.2 


0.13565 




0.4749(2) 


0.6794 (6) 


0.2045 (6) 


5.25 


0.1352 


0.54113(2) 








5.26 


0.1345 


0.53973(1) 


0.4739(4) 


0.6839(11) 


0.2100(11) 


5.29 


0.134 


0.54241(1) 


0.4707(3) 


0.6782(10) 


0.2075(10) 


5.29 


0.1350 


0.54552(3) 








5.29 


0.1355 


0.54708(3) 









Obviously many ways, some better than others, to ex- 
tract ctjjg exist that are consistent with perturbation 
theory at a given order. For convenience we adopt the 



Unfortunately, these results are of limited use since they have 
been obtained at re~ 1 values that correspond to negative quark 
masses, after subtracting re^T to the same order in perturbation 
theory. 



XVI 



P 



□ 



gym 



aV{2) 



F 12 



5.3 
5.3 
5.3 
5.3 
5.3 
5.3 
5.3 
5.35 
5.35 
5.35 
5.35 
5.35 
5.35 
5.35 
5.415 
5.415 
5.415 
5.415 
5.415 
5.415 
5.415 
5.5 
5.5 
5.5 
5.5 
5.6 
5.6 
5.6 
5.6 



0.3 
0.2 
0.15 
0.1 

0.075 
0.05 

0.025 
0.3 
0.2 
0.15 
0.1 

0.075 
0.05 

0.025 
0.3 
0.2 
0.15 
0.1 
0.05 

0.025 

0.0125 
0.1 
0.05 

0.025 

0.0125 
0.08 
0.04 
0.02 
0.01 



0.46980(6 
0.47554(8 
0.4792 (4 
0.48444(8 
0.48714(8 
0.48918(9 
0.49238(8 
0.48243(5 
0.48912(5 
0.49373(5 
0.49951(4 
0.50283(5 
0.50565(6 
0.51097(6 
0.50065(6 
0.50841(8 
0.51366(6 
0.52007(7 
0.52692(5 
0.53066(4 
0.5329 (1 
0.54328(2 
0.54724(2 
0.54932(2 
0.55027(2 
0.56325(2 
0.56418(1 
0.56479(1 
0.56499(1 



0.7055(35 
0.6848(34 
0.6740(17 
0.6534(24 
0.6512 (4 
0.6419(14 
0.6301(10 
0.6724(14 
0.6482(13 

0.6154(10 
0.6066(10 
0.5980 (9 
0.5822 (8 
0.6223(15 
0.5980(14 
0.5834 (8 
0.5627 (9 
0.5422 (5 
0.5295 (5 
0.5222 (5 
0.5070 (3 
0.4935 (1 
0.4862 (2 
0.4829 (2 
0.4616 (1 
0.4566 (1 
0.4537 (1 
0.4528 (1 



1.35 (14) 
1.114 (9) 
1.299 (6) 
1.291 (8) 
1.1580(18) 
1.135 (35) 
1.096 (20) 
1.258 (6) 
1.163 (36) 

1.024 (18) 
1.032 (18) 
1.014 (13) 
0.976 (10) 
1.113 (37) 
1.009 (15) 
0.977 (13) 
0.948 (8) 
0.8711(42) 
0.8449(35) 
0.8231(35) 
0.8001(17) 
0.7592 (5) 
0.7375(11) 
0.7289(11) 
0.6904 (7) 
0.6763 (5) 
0.6670 (7) 
0.6619 (6) 



0.64 (14) 
0.429 (9) 
0.625 (6) 
0.638 (8) 
0.5068(18) 
0.493 (35) 
0.466 (20) 
0.586 (6) 
0.515 (36) 

0.409 (18) 
0.425 (18) 
0.416 (13) 
0.394 (10) 
0.491 (37) 
0.411 (15) 
0.394 (13) 
0.385 (8) 
0.3289(42) 
0.3154(35) 
0.3009(35) 
0.2931(17) 
0.2657 (5) 
0.2513(11) 
0.2460(11) 
0.2288 (7) 
0.2197 (5) 
0.2133 (7) 
0.2091 (6) 



procedure detailed in Refs. [|[ but point out that 
many alternative ways are equally justified (e.g. that of 
Ref. ||). We only differ from Refs. [g, |j in as far as, once 
the coupling has been converted into the MS scheme, we 
evolve the scale and extract the A-parameter by numer- 
ically integrating the four-loop f3 function, rather than 
using a perturbatively truncated formula. 
We can write the plaquette as, 

□ = 1 - aa L - c 2 a 2 L - c 3 a 3 L , (89) 

where c\ and c 2 can be read off from Tab. [I] with the 
help of Eqs. (|l^) - (|||) for the pure gauge case and the 
fermionic contributions can be identified in Tabs. |IV| and 
In particular one finds c\ 



In particular one finds c\ = 47tCf/4. c 3 for the pure 
gauge case and sources in the fundamental representation 
reads M El: 



C3 



{N 2 - l)N 



0.0063538 



0.0181239 

N 2 



23 



0.0185221 



iV 4 



(4tt) 2 



We now define the quantity, 
lnD 



P 



Cl 



£2 
Cl 
£3 
Cl 



ci = — + — 



Co 



£l 

2 ' 



C2a|, 



(90) 



(91) 
(92) 
(93) 

(94) 



The expansion of P in terms of suffers from big coeffi- 
cients. For instance one obtains c\ ~ 3.37, c 2 ~ 17.69 in 
SU(3) pure gauge theory. The origin of these big num- 
bers can be traced to the lattice tadpole diagrams ||. 
Large coefficients are also encountered when aL is con- 
verted into the more "physical " cou pling ay [or otjjg 
which is "close" to ay, ci. Eq. (A55)]. By re-expressing 
the series Eq. ( p2[ ) in terms of ay, taken at a suitable 
scale q* jQ , one might hope that the two effects cancel 
in part |8|: 



P = a v (q*) + Xi(q*a)ay(q*) + X 2 (q*a)a v + ■ 



We obtain, 



X 1 {x) = h+B^x) 

X 2 {x) = c 2 + B 2 (x) + 2c- 1 B 1 (x), 



where 



B x {x) = 
B 2 {x) = 



—b\ — a\ + 2/3q In a; 

-b 2 - a 2 + 2(a 1 b 1 + b\ + a\) 

2/3i \nx + Bl{x)- B\{1). 



(95) 



(96) 
(97) 



(98) 
(99) 



The c onsta nt bi is defined in Eqs. (|A14|) - (|A17[) , b 2 in 
Eqs. (|A20|) - (|A2j) and ft in Eqs. QA2]) - (|A3|V The a; 
are defined in Eqs. ( A34) and (A37) and have to be taken 
at the scale q* for massive sea quarks. The numerical 
values for n f = and q*a ps 3.402 read: X 1 w —1.191, 
X 2 w -1.688. We can now truncate Eq. @ at C(a 2 ) 
and obtain, 



X 2 a 3 P 



2k 



1 - 1 



(100) 
(101) 



Note that the coefficients Xi inherit a mass dependence 
from a,-, &i (through AiQ) and c,-. With 



«ms(m) 


= a v (q*) - 


f *ia£(g*)- 




(102) 


/' 


= e-^ 6 q* 


ps 1.478/a, 




(103) 


ll 




r - 


§4 




(104) 


^2 




a 2 - 


«? - f A - 


(ai - §A>) 


C-05) 



we arrive at, 



where 



'MS 



Yx 
Y 2 



2N n f , 
In 

37T 67T 
2o — A"o. 



(106) 

(107) 
(108) 



For m = Yi is independent 9 of nj. For nf = 0, iV = 3 
we find the numerical value, Y2 ~ 0.9538: the NNLO 
correction is small. 

In analogy to defining a coupling from the measured 
plaquette, other couplings can be computed from force 
and potential. We illustrate this procedure at NLO: in a 
first step we can write, 



aV(Ra) 



aL 



u 2 (R) 
«i(R)' 



a v (q*) 



f 2 (R) 
fi(R) 



Bi{q*a) 



a 2 , (109) 



where the function B\(x) is defined in Eq. (pq). A 
similar expression can easily be written down for the 
force. The respective q*a values can be found in Tab. XI. 
Consequently, ay(q*) can be obtained by solving the 
quadratic equa tion E q. (|l Qj ) and converted into a-jj-g{p) 
via Eqs. (106) and (107), where we set Y 2 = in our 
0(a 2 ) calculation. 

Finally, we can run ajjg- numerically to arbitrarily high 
scales using the perturbative four-loop (3 function |63| , |64[j 
and then determine the Ajg^parameter, defined as in 



Eq. 



2. Pure gauge theory 



We display one- and two-loop results on ctjjg-(fi) as 
obtained from the logarithm of the average plaquette 
following the boosted perturbation theory proced ure de- 
tailed above (as well as in Refs. § 1), in Tab. |XXI1| 



This is then numerically converted into estimates of the 
QCD A^jg-parameter which are displayed in the last two 
columns of the table. In Fig. pB] we compare various 
methods of determining the M S A-parameter as a func- 
tion of the lattice spacing. The horizontal error band 
corresponds to the continuum limit result as obtained 
by the ALPHA Collaboration |n] by use of finite vol- 
ume techniques. aL refers to a conversion from the bare 
lattice coupling qj, = 3/[27r/3(a)] and ap to a coupling, 



9 For n f degenerate massive flavours one could in principle main- 
tain the mass independence of this coefficient: Yi 1— * Ki — 
llAT/(127r) ln(l + C m 2 /^ 2 ), fj, h-> fj, X yjl + C Q m 2 /fi 2 . In the 
most interesting case, QCD with non-degenerate flavours, this is 
however not possible. 
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"measured" from the average plaquette. NLO and NNLO 
refer to results from a naive perturbative conversion of 
aj, into ajgjfa" 1 ) at order a 2 and a 3 respectively while 
the abbreviations bNLO and bNNLO correspond to the 
boosted perturbation theory procedure detailed above. 



TABLE XXII: fj, = e~ 5/(i q* in SU(3) pure gauge theory and 
one- (1-1) and two- loop (2-1) estimates of o.-^-g{fi) and A^yg 
from the average plaquette. The errors of a do not reflect the 
uncertainty in the scale fi. 



P 



At 



,ro At 



5.5 
5.6 
5.7 
5.8 
5.9 
6.0 
6.2 



2.97 (4) 
3.61 (9) 
4.23 (7) 
5.38 (7) 
6.80(13) 
7.88 (4) 
10.78 (6) 

6.3 12.40(10) 

6.4 14.62(24) 
6.6 18.82(21) 



0.26368(3) 
0.22947(3) 
0.20417(3) 
0.18753(2) 
0.17582(2) 
0.16667(1) 
0.15233(1) 
0.14639(1) 
0.14105(1) 
0.13172(1) 



0.27517(3 
0.23747(3 
0.21000(3 
0.19216(2 
0.17970(2 
0.17000(1 
0.15494(1 
0.14873(1 
0.14316(1 
0.13346(1 



0.604 (9) 
0.568(15) 
0.517 (9) 
0.532 (7) 
0.566(12) 
0.562 (3) 
0.580 (3) 
0.584 (5) 
0.603(10) 
0.601 (7) 



0.648 (9) 
0.608(15) 
0.551 (9) 
0.566 (7) 
0.601(12) 
0.596 (3) 
0.613 (3) 
0.616 (5) 
0.636(10) 
0.632 (7) 



TABLE XXIII: Boosted NLO pure gauge estimates of cejrs(fi) 
from aV(l), aV(2), F12 and F23, respectively. Th e respe ctive 
scales (iro can easily be read off from Tabs. Kl and XVII . The 
errors of a do not reflect the uncertainty in [iro- 
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<4» 








5.5 


0.3414(17) 


0.629(27) 






5.6 


0.2862 (7) 


0.4409(32) 


1.212(39) 


1.594(89) 


5.7 


0.2443(13) 


0.3308(31) 


0.682(19) 


1.096(79) 


5.8 


0.2173 (6) 


0.2744(16) 


0.4958(90) 


0.721(30) 


5.9 


0.1984 (5) 


0.2426(10) 


0.4104(54) 


0.560(19) 


6.0 


0.1856 (1) 


0.2192(10) 


0.3444 (5) 


0.446(19) 


6.2 


0.1669 (1) 


0.1898 (5) 


0.2736 (3) 


0.332 (9) 


6.4 


0.1527 (1) 


0.1702 (2) 


0.2341 (8) 


0.271 (2) 


6.6 


0.1411 (1) 


0.1554 (1) 


0.2073 (6) 


0.226 (2) 



We find that the ratios between NLO and NNLO re- 
sults can be brought closer to unity, both by "measuring" 
the coupling and by "boosting" the perturbation theory. 
The combination of both methods indeed brings the re- 
sult in agreement with the known number (error band). 
Interestingly, almost all bNLO and bNNLO results are 
within the expected range, even at rather coarse lattice 
spacings. 



In Tab. XXIII we display values ajjg(fi) obtained from 
the short distance lattice potential and force 10 . These 
are then converted into estimates of A-parameters and 



TABLE XXIV: One-loop pure gauge estimates of A-jy^rrn from 
the average plaquette, aV(l), aV(2), F12 and F23, respec- 
tively. In the last line we display our continuum limit esti- 
mates. The errors are purely statistical. 



3 


A°-rn 


A Vl r n 






A F23 m 


5.5 


604 (Q) 


736(4 V[ 


QWCWl 






. D 


U.OOo^ 10 ) 




u.yoo^zy ) 






5.7 


0.517 (9) 


0.633(17) 


0.819(23) 








5.8 


0.532 (7) 


0.635(12) 


0.797(18) 


0.672(17) 


0.708(19) 


5.9 


0.566(12) 


0.653(18) 


0.813(23) 


0.734(24) 


0.806(32) 


6.0 


0.562 (3) 


0.639 (5) 


0.768 (5) 


0.708 (5) 


0.805 (8) 


6.2 


0.580 (3) 


0.650 (4) 


0.749 (4) 


0.703 (4) 


0.816 (7) 


6.3 


0.584 (5) 










6.4 


0.603(10) 


0.667(12) 


0.755(15) 


0.722(17) 


0.829(25) 


6.6 


0.601 (7) 


0.654 (9) 


0.736(10) 


0.717(13) 


0.764(24) 


00 


0.609 (4) 


0.666(5) 


0.727(8) 


0.735 (8) 


0.770 (9) 



compared to the corresponding bNLO estimates from 
the average plaquette in Tab. 



XXIII 



We also display 
these numbers versus the inverse momentum scale l//i 2 
in Fig. £26]. The curves represent results from phenomeno- 
logical quadratic plus quartic fits in l/(^tro). The 
Ajjg-ro values resulting from these continuum limit ex- 
trapolations are displayed in the last row of the table. 
The results from the high q* quantities V{\) and □ turn 
out to be in reasonable agreement with that from the AL- 
PHA Collaboration, A-p-gTn. = 0.602(48), however, this is 
not the case for Fn or V(2). 

Unfortunately, there exists no first principles way to 
estimate systematic uncertainties in any determination 
of Ajgij that is partially based on perturbation theory at 
momenta as low as a few GeV. For instance the F12 data 
exhibit a plateau for (3 > 5.9 that is in no way inferior to 
that obtained from the average plaquette. However, the 
two extrapolated values differ by almost three standard 
deviations of the ALPHA Collaboration result. This sug- 
gests that the good agreement obtained for the plaquette 
might be partly accidental. We are unaware of any con- 
vincing argument why an estimate of ajj^{fi) extracted 
from F12 at (3 — 6.6, i.e. {ir ~ 5.65 should be less reliable 
than the value obtained from the plaquette at say j3 = 5.8 
(/iro ~ 5.4). The same can be said about the combina- 
tion aV(2) at = 6.6, ^r « 12.8 and the plaquette at 



One could also imagine to repeat this procedure for a "tadpole 



improved" potential, V(R) +ln D/(2a) at the corresponding q*s. 
Since in this case the result is obviously correlated with that 
which we obtained from the plaquette we leave this exercise to the 
interested reader and remark that any such additive corrections 
cancel in the force. 

Since A cannot directly be obtained from position space Greens 
functions there is no theoretically well founded reason to assume 
the leading order lattice corrections to be quadratic in the lattice 
spacing. In fact we know that perturbative corrections to A exist 
which should be of order A/ ln[(aA) — 2 ]. 
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FIG. 25: Comparison of Ajg^g-parameters for SU(3) gauge 
theory. The error band denotes the continuum limit result ob- 
tained by the ALPHA Collaboration H . NNLO corresponds 
to 0{a 3 ), NLO to 0{a 2 ), bNLO and bNNLO to "boosted" 
NLO and NNLO, ctL to a conversion from the bare lattice cou- 
pling and ap to the a- values, "measured" from the logarithm 
of the average plaquette. 



observables as well as the ALPHA result at hand it be- 
comes obvious that the 4 % difference between bNLO and 
bNNLO estimates from the plaquette do not necessarily 
reflect the whole truth, but that a conservative approach 
would instead argue in favour of a 15 % error on the 
A-parameter in an 0(a 2 ) determination like the above. 
Quantities with higher q* values are certainly likely to 
behave better, but incorporating a number of measure- 
ments to estimate systematic errors is prudent when such 
a scale setting prescription has been used. 

The present article contains the NLO perturbation the- 
ory that is necessary for the use of additional short dis- 
tance quantities such as V(l) and the chair and par- 
allelogram Wilson loops. The observed scattering be- 
tween extrapolations from different input quantities will 
hopefully shrink both as finer lattices are simulated, and 
once higher order corrections are known for potential and 
force. In the meantime we quote the value, 

A^r = 0.609(4)(90), ( 110 ) 
where the second error is systematic. 



3. Results with sea quarks 
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FIG. 26: Comparison of Aj^g-parameters for 5(7(3) gauge 
theory from different non-perturbative input quantities as a 
function of the scale n = e~ 5 ^ 6 q* . The curves correspond to 
quadratic plus quartic fits in and the error band is the 
continuum limit ALPHA Collaboration result [fell. 



(3 = 6.3, firo w 12.4 and yet the former two (3 = 6.6 
estimates A|^Lr = 0.717(13) and A^r = 0.764(24) 
happen to lie significantly higher than the known contin- 
uum limit result. 

We conclude that without the comparison between pre- 
dictions from different non-perturbatively obtained ob- 
servables we might easily have underestimated the sys- 
tematic uncertainties of the approach. Having different 



We shall limit our discussion of the fcrmionic case to 
0(a 2 ) since only the plaquette is known to NNLO and 
only for Wilson quarks at some negative mass values [|l9j . 
The fermionic correction to C2 within Eqs. ( |89| ) - (p2| ) can 
be calculated from Tab. [v| with the help of Eqs. ( |l9| ) and 
the definitions Eqs. ( |l4| ) - ( fl6| ) . Th e corrections to the 
potential are available in Tabs. VIII - IX, with the con- 
ventions of Eqs. (p5|) and (pq) and v% 



(^yVi. We also 
include data that has been obtained by the CP-PACS 
Collaboration (3(J by combining SW fermions with the 
Iwasaki gluon action into the analysis. The correspond- 
ing perturbative expansion of the plaquette can be read 
off from Tabs. HI and IV, with the definitions Eqs. (LHj) 



- Q2ip. In this case h w -1.2466 + n f K? w (ma) |2J 
and K\ is independent of the gauge action. We find 
q*a w 3.213 for the plaquette. We display CP-PACS 
simulation para meter s and results on r^/a and the pla- 
quette in Tab. XXV . The quark masses are obtained 
via Eq. (|82j). We set csw = 1 in the analysis of SW- 
Wilson and SW-Iwasaki (SW-I) results which is consis- 
tent at NLO. 

We display our results on ajjg(fi) as well as on Ajjg-r$ 
in Tabs. IXXVll - IXXIXl. (5AE refers to the contribution 



to the A-p^-parameters due to the mass dependence of 
the coefficients b\ (through Ai\i) and is included into 
the Ajmj estimates. This mass dependence has not been 
considered prior to the present study. For SW and KS 
fermions the effect turns out to be of the order of the 
statistical uncertainty as long as ma < 0.1. However, for 
Wilson fermions this is a substantial effect, ranging from 
10 % at ma w 0.06 to 2 % at ma w 0.01. In this case the 
main effect is due to the term AK\(ma) in the matching 
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TABLE XXV: r /a and □ obtained on the CP-PACS ensem- 
ble |30|. 



P 


K 


ma 


ro/a 


1 — 1 
l_l 


1 on 
l.oU 


n 1 /t a a 
U.14U9 


a 1 fio 
U.162 


1. f lO^OO J 


n a one Q /q\ 
U.4y(Joo(oj 


1 on 
l.oU 


n 1 a q a 
U.14oU 


U.11U 


1. ryy^lo ) 


U.4yoUO(4j 


1.80 


0.1445 


0.074 


1 oh7 / a\ 

1.897(30) 


0.49936(4) 


l.oU 


U.1404 


n aoa 
U.U29 


2.U04(oo J 


U.OU ( 2U(0 ) 


1 A K 

1.90 


U.lo/0 


(1117 


2.49 I (04 J 


U.00oot>(2 ) 


1.90 


a 1 oaa 
U.139U 


U.O (8 


2.001(42 J 


0.0066/(2) 


1 a k 
1.90 


r\ t a a a 
U.14UU 


(1 nrn 

U.U02 


qoi ("OA 1 ! 
2.821(29) 


A £ K A1 /I /0\ 

U. 00914(2 ) 


1 a re 
1.90 


a 1 A 1 n 
U.141U 


n no77 
U.U2 < 


3.014(o3J 


0.06158(3) 


in 
2.1U 


a 1 9 c7 


a AQ7 


o.84o(lo) 


A KAQ AO / 1 \ 

U.09oUo(l) 


2.10 


0.1367 


0.060 


4.072(15) 


0.59920(1) 


2.10 


0.1374 


0.042 


4.236(14) 


0.60006(1) 


2.10 


0.1382 


0.020 


4.485(12) 


0.60108(1) 


2.20 


0.1351 


0.069 


4.913(21) 


0.62003(1) 


2.20 


0.1358 


0.050 


5.073(19) 


0.62062(1) 


2.20 


0.1363 


0.037 


5.237(22) 


0.62104(1) 


2.20 


0.1368 


0.023 


5.410(21) 


0.62149(1) 



TABLE XXVI: Boosted NLO estimates of o^sO) and 
AjjTg-ro from the average plaquette, with two flavours of Wil- 
son quarks. 5A^ and 5A^ J<1 refer to the corrections, included 
in Ajj-g, that are due to the finite quark masses and from the 
ma dependent term AKi alone, respectively. 



P 



5A° 



ro 



5A D 



5.5 0.060 

5.5 0.040 

5.5 0.028 

5.5 0.020 

5.6 0.050 
5.6 0.040 
5.6 0.030 
5.6 0.020 
5.6 0.010 



5.95 (4) 
6.48 (4) 
6.91 (5) 
7.23 (4) 
7.54 (4) 
7.81 (8) 
8.09(11) 
8.71 (4) 
9.21 (9) 



0.20925(1 
0.20475(2 
0.20103(3 
0.19988(5 
0.19405(1 
0.19247(1 
0.19076(1 
0.18885(1 
0.18702(1 



0.549(3) 
0.564(4) 
0.578(4) 
0.588(4) 
0.564(3) 
0.570(6) 
0.575(8) 
0.601(3) 
0.617(6) 



+0.049 
+0.038 
+0.030 
+0.023 
+0.045 
+0.038 
+0.031 
+0.023 
+0.012 



+0.056 
+0.043 
+0.034 
+0.026 
+0.051 
+0.043 
+0.035 
+0.026 
+0.013 



between lattice and MS scheme and we hence conclude 
that un-improved Wilson fermions are not particularly 
well suited for unambiguous determinations of the run- 
ning coupling, unless data at rather small quark masses 
are available. We have neglected the mass d ependence 
in our analysis of the SW-I data (Tab. XXIX), which is 
equally small as in the SW- Wilson case. 

We use our perturbative result on the potential to es- 



timate the ratio Ajj-gTo from the potential aV(l), aV(2) 
and force F12 for the three rif = 2 data sets with Wilson 
glue. The resulting numbers are compared t o the results 
from the plaquette in Tabs. [XXX| - |XXXIl| . Like in the 
quenched case the estimates from V(2) consistently turn 
out to be bigger by about 15 %, in comparison to those 
obtained from 1^(1) or the plaquette, however, in the case 



TABLE XXVII: Estimates of «i7s(m) an d JAws r o from the 
average plaquette, with two flavours of SW quarks. SA^ refers 
to the finite quark mass correction, included in A^jg. 






ma 








8A°r 


5.2 


0.046 


7.03 (6) 


0.19109(1) 


0.502(4) 


+0.001 


5.2 


0.024 


7.45 (6) 


0.18903(1) 


0.515(4) 


+0.004 


5.25 


0.043 


7.59 (7) 


0.18568(1) 


0.497(5) 


+0.005 


5.26 


0.072 


6.96 (8) 


0.18647(3) 


0.462(5) 


+0.003 


5.29 


0.093 


7.11 (7) 


0.18407(5) 


0.454(5) 


+0.000 


5.29 


0.053 


7.78(10) 


0.18255(1) 


0.483(7) 


+0.005 


5.29 


0.035 


8.31(13) 


0.18143(1) 


0.507(8) 


+0.005 
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FIG. 27: Comparison between bNLO ajfs estimates from 
pure gauge theory simulations and simulations with nj = 2 
Wilson, KS and SW quarks, the latter with Wilson and 
Iwasaki (SW-I) gluon actions. The n/ = error band cor- 
responds to the ALPHA result A w = 0.6Q2(48)r ( T 1 while 
the rif = 2 band corresponds to the four- loop running with 
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0.69(15). 



of KS quarks the estimates from Fyi are somewhat more 
in line with those from the plaquette. At present we do 
not have data at sufficiently many different lattice spac- 
ings for a detailed comparison like the one presented in 
Fig. ^6] for the n/ = case. Hence, we will only use 
the most ultra violet quantity, the average plaquette, to 
extract the nf = 2 A-parameter and assume a system- 
atic uncertainty of 15 % due to higher order perturbative 
corrections, based on our study of the p ure gaug e case as 
well as on the numbers in Tabs. XXX - XXXII. 

In Fig. we compare our nf — 2 estimates on ajjg(fi) 
from the plaquette to the nj = results. Due to the vari- 
ation in the sea quark masses, action and lattice spacings 
used, the nf = 2 points scatter around quite a bit but 
yield consistently larger values of a, compared to the 
quenched case as already observed in lattice simulations 
of the static potential, for instance in Ref. |39t] as well as 
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TABLE XXVIII: Estimates of a w (/i) 
average plaquette, with "two" flavours 



and A-jy^ro from the 
of KS quarks. 



TABLE XXIX: Estimates of cejj-g(fi) and A-j^rro from the av- 
erage plaquette, with two flavours of SW quarks and Iwasaki 
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ro 
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u.ou 
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Z. 44(^4 J 


U.OODZ (ioj 


0.587(15) 


-u.uoz 


E QE 

O.oO 


n qa 
U.oU 


O fiE/o\ 
Z.OO^Z j 


A Q A£Q (1 A \ 
U.oUoy [14J 


0.552 (7) 
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u.ou 


O A1 
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0.496(16) 
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0.19960 (3) 


0.575(12) 


-0.001 


5.415 


0.0125 


4.54(3) 


0.23492(13) 


0.558 


(4) 


-0.000 


5.5 


0.0125 


5.88(2) 


0.21447 (2) 


0.579 


(2) 


-0.000 


•3.6 


0.08 


6.03(2) 


0.20033 (3) 


0.494 


(2) 


-0.004 


5.6 


0.04 


6.71(3) 


0.19975 (1) 


0.545 


(3) 


-0.002 


•3.6 


0.02 


7.15(2) 


0.19929 (1) 


0.577 


(2) 


-0.001 


5.6 


0.01 


7.38(2) 


0.19916 (1) 


0.594 


(2) 


-0.000 



in Sec. IV C| below. It is also clear from the figure that at 
present the /i- window covered by rif — 2 results is much 
smaller than the one available in pure gauge studies. 

Since chiral symmetry is explicitly broken for the 
fermionic actions used we should first extrapolate our 

(2) 

Aij^. estimates obtained at similar quark masses in phys- 
ical units mro to the continuum limit, before attempting 
a chiral extrapolation: m — > m Ul m e i « 0. Due to the 
small range of lattice spacings covered by the presently 
available data, however, we are forced to interchange the 
ordering of these two limits: after discarding lattices with 
2.5 a > ro, i.e. the KS data with < 5.415 or ma > 0.05, 

(2) 

we extrapolate in Aj^rg at fixed p- values to small quark 
masses and subsequently send the lattice spacing to zero. 

The chiral extrapolations arc illustrated in Fig. ^ for 
Wilson glue and in Fig. ^9] for the SW-I action. We in- 
clude the /3 = 1.8 data points into the latter figure but we 
will discard these coarse lattice results from our contin- 
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2.10 


0.042 


5.915(20) 


0.21790(1) 


0.606 (2) 


2.10 


0.020 


6.263(17) 


0.21729(1) 


0.637 (2) 


2.20 


0.069 


6.860(29) 


0.20616(1) 


0.608 (3) 


2.20 


0.050 


7.084(27) 


0.20582(1) 


0.625 (3) 


2.20 


0.037 


7.313(31) 


0.20558(1) 


0.643 (3) 


2.20 


0.023 


7.554(29) 


0.20531(1) 


0.662 (3) 



TABLE XXX: Boosted NLO estimates of A^r for Wilson 

M S 

quarks from the average plaquette, aV(l), aV(2) and -F12, 
respectively. The errors are purely statistical. 



p 


ma 


A w r ° 


4f^o 


^0 




5.5 


0.060 


0.549(3) 


0.616 (5) 


0.786 (8) 


0.687 (8) 


5.5 


0.040 


0.564(4) 


0.628 (5) 


0.784 (8) 


0.690 (8) 


5.5 


0.028 


0.578(4) 


0.641 (6) 


0.790 (8) 


0.697 (8) 


5.5 


0.020 


0.588(4) 


0.650 (6) 


0.788 (8) 


0.693 (8) 


5.6 


0.050 


0.564(3) 


0.621 (5) 


0.759 (7) 


0.691 (7) 


5.6 


0.040 


0.570(6) 


0.627 (8) 


0.758(11) 


0.686(11) 


5.6 


0.030 


0.575(8) 


0.634(10) 


0.755(14) 


0.678(15) 


5.6 


0.020 


0.601(3) 


0.659 (4) 


0.784 (6) 


0.709 (7) 


5.6 


0.010 


0.617(6) 


0.675 (8) 


0.805(13) 


0.734(15) 



uum limit extrapolation. The last column of Tab. XXX 
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reveals that the slopes of AjJ^ro as a function of mro 
would have been much steeper for Wilson fermions had 
we neglected the ma dependence of the matching be- 
tween lattice and MS schemes. In the course of each 
extrapolation the lattice spacing in physical units a/ro 
changes. However, A is no spectral quantity and there is 
no theoretical handle on the functional form of such an 
extrapolation anyway. We find linear fits, 



4^ 



{mro) 



(0) + cmro, 



(111) 



(2) 

with parameters [A-L_Lro](0) and c, to be consistent with 
the mr < 0.25 (mr < 0.35 for the SW-I action) data. 
The resulting parameter values as well as a/ro(m = 0) 
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TABLE XXXI: Estimates of AgLr for SW quarks. 
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ma 


* n 

A W r ° 


a Vi 


A if^o 




5.2 


0.046 


0.502(4) 


0.537(5) 


0.653 (8) 


0.672 (9) 


5.2 


0.024 


0.515(4) 


0.546(5) 


0.655 (7) 


0.670(10) 


5.2 


0.018 




0.559(6) 


0.666 (8) 


0.678(10) 


5.26 


0.072 


0.462(5) 


0.500(7) 


0.611(10) 


0.645(12) 


5.29 


0.093 


0.454(5) 


0.496(6) 


0.605 (8) 


0.648(11) 



TABLE XXXII: Estimates of A^r for KS quarks. 



p 


ma 


AS_ 


ro 


. M a . 


A w r ° 


A ^o 

. ad a 


5.3 


0.300 


0.587(15) 


0.684(19) 






5.3 


0.200 


0.603(19) 


0.699(27) 






5.3 


0.150 


0.595 


(5) 


0.694 (8) 






5.3 


0.100 


0.592 


(6) 


0.683(13) 






5.3 


0.075 


0.592 


(7) 


0.699 (8) 






5.3 


0.050 


0.596 


(7) 


0.696(11) 






5.3 


0.025 


0.697 


(4) 


0.692 (6) 






5.35 


0.300 


0.552 


(7) 


0.670 (7) 






5.35 


0.200 


0.551 


(5) 


0.661 (7) 






5.35 


0.100 


0.546 


(3) 


0.650 (6) 


0.840(31) 




5.35 


0.075 


0.556 


(3) 


0.664 (6) 


0.900(33) 




5.35 


0.050 


0.590 


(8) 


0.702(12) 


0.959(34) 


0.434(71) 


5.35 


0.025 


0.559 


(5) 


0.660 (8) 


0.904(25) 


0.526 (9) 


5.415 


0.300 


0.496( 


16) 


0.618(19) 


0.879(76) 




5.415 


0.200 


0.508 


(4) 


0.626 (7) 


0.846(26) 


0.477(13) 


5.415 


0.150 


0.510 


(3) 


0.627 (5) 


0.854(25) 


0.517(12) 


5.415 


0.100 


0.515 


(2) 


0.620 (5) 


0.881(18) 


0.561 (7) 


5.415 


0.050 


0.530 


(2) 


0.625 (4) 


0.836(11) 


0.596 (7) 


5.415 


0.025 


0.551 


(2) 


0.634 (4) 


0.853(11) 


0.630 (7) 


5.415 


0.0125 


0.558 


(4) 


0.634 (6) 


0.840(13) 


0.635(10) 


5.5 


0.100 


0.488 


( 7 ) 


0.580 (9) 


0.787(14) 


0.623(11) 


5.5 


0.050 


0.515 


(5) 


0.592 (6) 


0.772 (8) 


0.631 (7) 


5.5 


0.025 


0.559 


(3) 


0.629 (4) 


0.800 (7) 


0.661 (7) 


5.5 


0.0125 


0.579 


(2) 


0.645 (3) 


0.814 (5) 


0.676 (5) 


5.6 


0.080 


0.494 


(2) 


0.564 (2) 


0.717 (4) 


0.627 (4) 


5.6 


0.040 


0.545 


(3) 


0.607 (3) 


0.756 (5) 


0.662 (5) 


5.6 


0.020 


0.577 


(2) 


0.634 (2) 


0.776 (4) 


0.676 (4) 


5.6 


0.010 


0.594 


(2) 


0.650 (2) 


0.783 (4) 


0.676 (4) 



TABLE XXXIII: Chirally extrapolated a/ro values and ex- 
trapolations of A-^L-ro, according to Eq. (111). 



Action 


p 


a/ro 


c 


[ A w^](°) 


x/Ndf 


Wilson 


5.5 


0.178(4) 


-0.271(17) 


0.613(3) 


0.51/2 


Wilson 


5.6 


0.148(3) 


-0.335(47) 


0.639(7) 


1.67/2 


SW 


5.2 


0.185 


-0.136 


0.532 


Ndf = 


SW 


5.29 


0.162(4) 


-0.198(23) 


0.542(9) 


0.48/1 


KS 


5.415 


0.312(2) 


-0.311(30) 


0.573(3) 


0.56/1 


KS 


5.5 


0.236(6) 


-0.508(38) 


0.605(3) 


1.06/1 


KS 


5.6 


0.194(1) 


-0.371 (8) 


0.613(1) 


0.16/1 


SW-I 


1.8 


0.460(11) 


-0.251(34) 


0.508(7) 


0.78/2 


SW-I 


1.95 


0.312 (5) 


-0.358(34) 


0.575(6) 


0.83/2 


SW-I 


2.1 


0.213 (1) 


-0.327(12) 


0.666(3) 


2.08/2 


SW-I 


2.2 


0.177 (1) 


-0.254(18) 


0.692(4) 


2.32/2 
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FIG. 28: Chiral extrapolations of Ai2^.ro at finite lattice spac- 
ings for the Wilson gluonic action combined with the three 
different fermionic actions. 



0.75/1 for KS quarks, A^r 



w- - 0J64 and A ws r ° - 

0.575 for Wilson and SW fermions, respectively, in which 

,(2) 



0. For the SW-I action we obtain Ai-2_ r n 

MS u 

0.747(2) with x 2 /N D f = 0.05/1 



cases N DF 

As a final result we 
quote the (unweighted) average of the KS and SW-I data 
sets, 



can be found in Tab. pCXXIIl| . For the SW data at 

and at f3 
0.45. 



5.29 we 



5.20 no fit is possible since Ndf 
included a rather heavy mass, mro 

Finally we attempt a continuum limit extrapolation 
in Fig. ^0|, where we have generously estimated error 
bars for the SW quarks at j3 = 5.20. The lines repre- 
sent extrapolations, linear in a/ro f° r Wilson fermions 
and quadratic in a/ro f° r KS and SW quarks. The ex- 
trapolated values are A^Lr = 0.638(3) with x 2 ' /Ndf — 



A Ws r ° = 0.69 ±0.15, 



(112) 



where the error estimate reflects the uncertainties due 
to higher order perturbative corrections as well as in the 
extrapolations, both chiral and to the continuum limit. 

The above number should be related to the recent more 
optimistic estimate A^Lr = 0.553(34) by the QCDSF 
and UKQCD collaborations M which is based on the 
UKQCD SW data alone. Note that their result agrees 

(2) 

with our above value, Al^ro ~ 0.575, obtained from the 

' MS v ' 
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FIG. 29: Chiral extrapolations of Ai^Lro at finite lattice spac- 
ings for the SW-I action. 
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FIG. 30: Attempts of a continuum limit extrapolation (linear 
in a/ro for Wilson and quadratic for KS and SW quarks) of 
A— — rn from the available "world lattice data" . 

MS 



same SW data set by use of somewhat different methods 
in the conversion between schemes and extrapolations. 
For instance in Ref. || cejjg(fi) was extracted at NNLO, 
estimating the as yet unknown NNLO fermionic contri- 
bution to the plaquette, while we consistently worked at 
NLO. 

Fig. |3(] paints a superficially disappointing picture for 
universality of the continuum limit. Universality abso- 
lutely must hold for lattice field theory to be a worth- 
while exercise, so one might be concerned about the large 
difference in particular between the two different results 
obtained with SW fermions. While we try to estimate 
systematics by incorporating this difference into our er- 
ror, it is interesting to consider how one could expect the 



figure to mature as improved data from future simula- 
tions become available. 

Firstly, we note that the KS and SW-I continuum limit 
extrapolations are much better controlled than those 
for Wilson or SW quarks with Wilson glue, due to the 
availability of an additional data point. There are siz- 
able 0(a 2 ) corrections to the current KS and SW data. 
The Wilson continuum limit result has been obtained 
by drawing a straight line through two data points but 
there is no reason to believe 0(a 2 ) corrections to be any 
smaller (nor indeed any larger) in this case than in the 
SW and KS cases. Likewise, it is not clear whether part 
of the large slope of the SW-I data set (which is only 
perturbatively improved such that 0(ar,a) corrections 
are still present) might be due to a residual linear con- 
tribution. One might hope that the scatter of the results 
will be reduced by future data at different lattice spac- 
ings that would allow us to discriminate between linear 
and quadratic terms. 

Secondly, we note that perturbative coefficients for the 
effects of the clover term in the SW action are particu- 
larly la rge th roughout the literature, and refer the reader 
to Eq. ( A16 ) for an ex ample of an 0(a 2 ) coefficient and 
Eqs. (|A23|) and ( |A24j ) for 0(a 3 ). There is no automa- 
tism that guarantees 0(a 3 ) corrections to the plaquette 
in determinations of otjjg to be small, just because this 
happened to be the case in the quenched theory with 
Wilson glue. Given the experience of the large c v sw cor- 
rections we find it at least plausible that perturbative 
corrections to the SW determinations are larger than for 
the two other formulations. A calculation of these might 
result in a better agreement between the different con- 
tinuum extrapolations too. 

We remain confident that a combination of more lat- 
tice data and higher order perturbative calculations will 
greatly reduce the systematic spread shown in Fig. ||(]. 
Of course eventually the ordering of chiral and continuum 
limit extrapolations should be interchanged, unless data 
with chirally symmetric lattice fermions become avail- 
able. 

We finally note that efforts are on the way of calculat- 
ing the rif — 2 A-parameter with the help of finite size 
scaling techniques in the Schrodinger functional scheme. 
At present the error from this approach |6(| is of O(50%) 
and the scale has still to be related to a physical quan- 
tity like ro or however, some improvement can be 
expected in the near future. 



4- Implications for the "real" world 

(2) 

Running the estimate A-^-ro ~ 0.69(15) to the mass of 
the Z vector boson mz ~ 91.2 GeV results in the value, 



o|( m z) = 0.0945(30)(7), 



(113) 



where the second error that is smaller than 1 % is due 
to the 5 % scale ambiguity |67l| in the phenomenological 
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value ro ~ 0.5 frn. This subdominant uncertainty can be 
reduced by using a quantity directly accessible in exper- 
iment like /„■ or mjv to set the scale. Nature of course 
does not contain two but five quark flavours lighter than 
the Z boson and we have to address this problem if we 
intend to produce numbers that have phenomenological 
relevance. 



0.13 
0.125 
0.12 



N 

E. 

£ 0.115 

CO 

a 

0.11 



0.105 



Lattice ' — ■ — i 
PDG ^—-^ 
Eq.(115) 




~n — 






i 











0.1 

[1] [2] [3] [4] [5] [6] [68] 
Reference 

FIG. 31: Our estimate of a^|L (error band), in comparison 
with numbers quoted in previous lattice studies as well as with 
the average from the "Review of Particle Physics" (PDG). 



The main difficulty of a QCD prediction is the effect 
of the strange quark. Given that we cannot resolve the 
shift between A^Lr = 0.60(5) and A^r = 0.69(15) 
within the errors, for the moment being we assume that 
the strange quark which is heavier than the up and the 
down quarks anyway will not have a large impact on this 

(3) (2) 

ratio either and guess A— — rn w AXJ—rn- However, we 

° MS MS 

try to incorporate the uncertainty of this assumption (as 
well as the scale error of ro) by somewhat inflating the 
systematic error: 



= 270(70) MeV. 



(114) 



Our result of Sec. 
tive matching of the 



[II E indicates that a perturba- 
rif — 3 and nt 



4 couplings 



around the charm threshold should be reliable. Fur- 
ther, the difference between a three and a four- loop MS 
running from a scale of 1.1 GeV upwards corresponds 
to a relative shift in a-jjg at the Z resonance of only 
0.26 %. Hence, we perturbatively match the coupling at 
the charm and bottom thresholds. We allow pole masses 
1.1 GeV< m c < 1.5 GeV and 4.4 GeV < m b < 5.0 GeV 
and arrive at the estimate, 



between a three and four-loop running of the coupling 
have been added linearly in the above error. 

The result is compatible with the average quoted in 
the "Review of Particle Physics" p8| : a^(mz) = 
0.1172(20) but lower than results from DO/CDF or LEP 
experiments, a tendency that is consistent with previ- 
ous lattice estimates 0, |, || J§ §. In Fig. |l| we re- 
late various lattice studies in chronological order to ours. 
We find that, although all lattice results on this quan- 
tity JE |, |, § |, § (as well as the PDG world aver- 
age [p8| ) lie within the error bars that we obtain, many of 
the previous numbers turn out to be incompatible with 
each other, within their quoted uncertainties. For in- 
stance we obtain a y^/NoF = 17.2/5 when averaging 
the six lattice results, and this assuming that the quoted 
errors, that are dominated by systematics, only represent 
one standard deviation! 

The least well controlled step in our prediction 

Eq. (|ll5| ) is certainly the extrapolation n/ = 0,2 — > 
nf = 3. A first lattice study with "n/ = 3" exists [||, 
however, it is too early to investigate quark mass effects 
and to attempt a continuum limit extrapolation in this 
case. Once lattice results with 2+1 quark flavours of a 
similar quality as those analysed here become available 
a determination of ctjfg(mz) with about 4 % accuracy 
using the methods explored above should become real- 
istic, and by both, a full fledged 0(a 3 ) calculation in 
lattice perturbation theory and more simulation points 
that would allow for a better controlled continuum limit 
extrapolation this error should be reducible to 1 - 2 % 
within the next few years. 

We experienced that systematic errors are easily un- 
derestimated in this kind of study, in particular since 
neither the size of higher order perturbative corrections 
nor the functional form of the chiral and continuum limit 
extrapolations are fully controlled. However, using a set 
of different lattice observables and actions certainly helps 
in arriving at realistic error estimates. Perturbation the- 
ory is naturally most reliable when the scales involved are 
high and it is very expensive to increase the energy a -1 
in lattice simulations, at least if one keeps the simulation 
volume fixed in physical units. This makes the average 
plaquette a particularly valuable quantity in studies like 
the present one. Future predictions might benefit if lat- 
tice data on the plaquette, not only in the fundamental 
but also in the adjoint and higher representations, were 
available as well as data on other short distance quanti- 
ties like the "chair", the "parallelogram" and the "rect- 
angle" . 



Comparison to the non-perturbative static 
potential 



( 5 ) 1 \ 
a MS^ 



0.1133(59), 



(115) 



using the four-loop /3-function. The effect of varying the 
respective quark masses m c and as well as differences 



We shall compare our perturbation theory results with 
data from simulations of QCD without and with sea 
quarks. First we explain what we mean by boosted 
and tadpole improved boosted perturbation theory in 
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the present context. Next we investigate different pos- 
sibilities to parametrize the pure gauge static potential 
at short distances, based on our results of Sec. III. We 
conclude with an attempt to resolve the different run- 
ning of the coupling when including sea quarks from non- 
perturbative data. 



Parametrizations of the potential 



In what follows we will refer to the parametrization 



aV(Ra) = ff 2 Vt(R) 



(116) 



with V\ defined in Eq. ( |30| ) [and Eq. ( p3 )] as leading order 
(LO). Apart from a different value of Ijr = 6//3 this result 
does not depend on nj. The NLO result can be found in 
Eqs. (|3l| ) - (|3^). The LO result above corresponds to 



aV(Ra) - aV s (a) 



-C* 



AttR 



(R 



(117) 



where aVs = CfQ 1 x 0.252731 . . . while the parametriza- 
tions of the NLO and NNLO curves in the R 3> 1 limit 
with an d w ithout massive sea quarks can be found in Ap- 
pendix A_4. In the pure gauge case we will also compare 
against NNLO expectations for which the small R lattice 
effects have not yet been calculated. When including 
fermions this is not possible since even in the massless 
case the contribution to the self energy is only known to 
NLO. 

In addition to expansions in terms of the bare cou- 
pling we shall also compare against "boosted" perturba- 
tion theory exp ectatio ns. In this case we use the A 
values of Tabs. |XXIV 



and 



XXVI 



MS 

calculated from the 
Since we shall 



respective plaquette at NLO as input 
compute the coupling ay(q) at NNLO this is somewhat 
inconsistent. However, we also aim at a comparison with 
the fermionic case where such a determination of the A- 
parameter is at present only possible at NLO and prefer 
to compare like with like. As we have seen above, the 
numerical difference between NLO and NNLO turns out 
to be small anyway for this quantity. Using A^g^ro as 
input we numerically run ctjfs(n) to the required scales 



/i = e~ 5 / 6 g* , according to the four-loop /3-function. Some 
of the respective aq*(R) values can be read off Tab. XI 
and for large R we employ t he par ametrization Eq. (60). 
Slightly deviating from Sec. [V B , the resulting ajjg(p) 
is then always converted into av(q*) at the two loop 
level [<D(g*% at bLO as well as at bNLO and bNNLO, 
via Eqs. flA33| ) - ( |A43[ ) and Eq. ( |A47[ ). This means that 
what we call bLO includes some NLO and NNLO contri- 
butions and bNLO includes a conversion from the MS to 
the V scheme at NNLO. However, the conversion from 
the momentum space continuum to the position space 
lattice potential is done at the order indicated. 

Finally, we introduce "tadpole improved" boosted per- 
turbation theory. We start from the observation that 



each link that appears within any lattice observable con- 
tains classes of tadpole diagrams that are lattice spe- 
cific. We try to cancel these ultra violet contributions in 
part by a non-perturbative prescription (?], B]: the sec- 
ond of the two Wilson loops that appear within the dif- 
ference aV(aR) = lim T ^ 00 {ln[M/(R,T)] - ln[W(R,T + 
1)]} contains two more links than the first one. We 
hence prefer to perturbatively expand the combination 
aV(aR) + lnD/2 instead, within which the net number 
of links is zero. We have already made use of this trick 
in Eq. (|74|). Finally, we reexpress the result in terms of 
cky(g^r), boosting our perturbative series, where the new 
(R) values can easily be calculated from those of the 
potential q* (R) and the plaquette q^ : 



Haq* N (R)) = 



ln[a g *(R)]Fi(R)- 111(^5)01/2 
Vi(R)- Cl /2 



(118) 



The av(q%) are then obtained as described above. We 
shall refer to this third method as tbLO, tbNLO and 
tbNNLO, depending on the order of the conversion from 
a v [q* N (R)} to aV(R). 



2. Pure gauge theory 

In Fig. H we compare LO, NLO, NNLO and bLO, 
bNLO and bNNLO expectations against the static po- 
tential at j3 = 6.0. At this lattice spacing we use 
Ajfjro = 0.562(3) and there is no free parameter. The 
un-boosted results come closer to the non-perturbative 
data as the order of the calculation increases but even at 
NNLO only about half of the potential can be explained, 
even at distances as short as R = 1. The boosted re- 
sults are more in line with the data but still only on a 
very qualitative level at best. It is also obvious that if 
one tuned the A-parameter such that V(l) is reproduced, 
V(2) would be undershot, in agreement with our expe- 
rience from Sec. IV B above. We can define an effective 



Coulomb coupling, 



e(R) = -R[aV(R) - aV s ] 



(119) 



and find e w 0.106 at LO while e(3) « 0.19 and e(3) « 
0.26 at NLO and NNLO, respectively. Finally, bNLO 
results in e(3) « 0.49. 

Next in Fig. || we try out the tbLO, tbNLO and 
tbNNLO parametrizations and indeed tadpole improve- 
ment vastly reduces the difference with the NP potential. 
To exclude that this is just accidental we also display a 
comparison with the potential at j3 = 6.4 in Fig. |^ where 
the plaquette yields Ajj^tq — 0.603(10), and with the 
smaller gauge coupling the situation becomes even more 
convincing. It turns out to be absolutely essential to com- 
pare like with like. Contrary to previous claims p0| , |69| ] 
not only the force but also the short range potential can 
be understood in terms of perturbation theory, provided 
that the self energy contribution is dealt with in a con- 
sistent way: when comparing with NP lattice results the 
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FIG. 32: Comparison between quenched lattice data at 
P = 6.0 and boosted (b) and un-boosted LO, NLO and NNLO 
perturbative potentials. The curves are the large 7? expecta- 
tions from continuum perturbation theory. 



FIG. 34: The same as Fig. B3l at (3 = 6.4. 



was found. 



CCS 

QC 
> 

CO 



1 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 



comparison with boosted and tadpole improved perturbation. 



self energy should not be subtracted from the perturba- 
tive expansion while when matching with the MS scheme 
in which Vg* s = by definition it has to be subtracted 
at the same scales fi(r) at which the interaction energy 
Vi n t(r) = V(fi;r) — Vs{p) is calculated, to avoid renor- 
malon ambiguities. A similar observation has been made 
in a recent comparison with continuum perturbation the- 
ory in a renormalon based approach [frof . In Ref. |^9[ a 
comparison was made between the NP lattice potential 
at short distances and the perturbative MS scheme pre- 
diction. In this study one and the same perturbative self 
energy had been used at all distances, which was obtained 
from a fit to the data, and consequently, no agreement 
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FIG. 35: Non-Perturbative minus pertur bativ e potentials at 
/3 = 6.4, together with a quadratic fit, Eq. ( |l20| ) . The distance 
r = 0.4 ro « 1 GeV - corresponds to about 4 lattice spacings. 



In Fig. |35| we display the difference between non- 
perturbative and perturbative potentials, AV(r) — 
V(r) — T/tbNLO^-j £ or ^.jjg exam pi e of /3 = 6.4 as a func- 
tion of R — r/a in physical units r^/a = 9.89(16). Within 
the statistical errors of the non-perturbative potential all 
points lie on a smooth curve indicating that the lattice 
effects are accounted for by the perturbative result. How- 
ever, there arc deviations between prediction and lat- 
tice data: at short distances we overestimate the NP 
potential by 30 MeV and around 0.3 fm we underesti- 
mate the result by 110 MeV. We fit the difference for 
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r < 0.61 ro w 0.3 fm by the phenomenological ansatz, 

AV(r) = b + cr 2 (120) 

with b = -0.068(3) Tq 1 , c = 0.99(6) r^ 3 and x 2 /Ndf = 
18.2/17. Allowing for an additional linear term does not 
improve the quality of the fit and the term turns out to be 
in agreement with zero: unlike in Ref. |69 no linear term 
is required at short distances to understand the data and 
this without any free parameter in the computation of 
AV(r)\ We notice that the coefficient c w (1.6A^jg) 3 is 
of order one in units of the QCD A-parameter. 




FIG. 36: The running of the Coulomb coupling at /3 = 6.6, 
where 5.1a « 0.4r « 0.2 fm « 1 GeV -1 . We have sub- 
tracted the value aVg hN w 0.65 from the non-perturbative 
potential. 
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FIG. 37: The same as Fig. || at /3 = 6.0, where 2.1 a « 0.4 
and we subtract the estimate aVg ~ 0.83. 
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In view of our plan to resolve the differences in the run- 
ning of the coupling between quenched and un-quenched 



simulations, in Fig. p6| w e dis play the dimensionless com- 
bination e(r/a) of Eq. flll9| ) for (3 — 6.6 where 5.1a ps 
0.4 r ps 1 GeV" 1 . At this /? value the average plaquette 
yields A-p^ro = 0.601(7). We have also included bare 
LO, NLO and NNLO results for comparison. The self 
energies have been subtracted from all perturbative po- 
tentials, LO, NLO, NNLO as well as tbNLO at the re- 
spective orders in perturbation theory. As what we are 
doing is a linear transformation we do not wish to sub- 
tract different V$ values from the NP and tbNLO poten- 
tials which would spoil the qualitative agreement between 
the results found above. Ideally one would subtract a 
non-perturbatively determined Vs which then guarantees 
scaling between different NP data sets. The difference 
between NP and tbNLO however will then result in the 
tbNLO coupling to diverge linearly at large R where per- 
turbation theory will become unreliable anyway. In the 
absence of a non-perturbative determination of Vs we will 
take the f3 = 6.6 tbNLO value aV| bNLO « 0.650 as an 
estimate which is indeed close at least to the tbNNLO 
result aVj bNNLO ps 0.670. Doing this means that the 
tbNLO Coulomb coupling will increase logarithmically 
with R and no linear term is present. 

In the figure we see a gap opening up between NP 
and tbNLO potentials for Ra > 4a ps 0.15 fm. The NP 
data starts diverging towards negative values, which is 
expected as e(r) — * — or 2 as r — > oo where a denotes 
the string tension. In addition, our estimate of an NP 
self energy might be wrong by a few per cent which will 
result in an unwanted linear contribution to e. We find 
relatively large values of the Coulomb coupling, e ps 0.4, 
that naturally disagree with e ps 0.3 as obtained from 
phenomenological fits to a Cornell type potential Eo| , 
aV(R) = Vo — e/R + a 2 a R, in which the linear confin- 
ing term contributes at all distances: the parameter e 
within this parametrization can only be interpreted as 
a Coulomb coupling at distances r <C \J~5. We observe 
that the running of the coupling is quite steep at short 
distances which provides us with hope that differences 
between the quenched and un-quenched /3-functions can 
be resolved from simulation data. 

In Fig. |37] we display the same comparison for a 
similar R range in physical units at (3 = 6.0 where 
2.1a w 0.4 r ~ lGeV" 1 . In this case the differ- 
ence between tbNLO and tbNNLO self energies is 15 % 
such that we cannot regard these values as reliable any- 
more. Instead we compute a 60 Vg ° — a 60 V 6 -°(ro) — 
(R%- 6 /R%-°)a 6 - 6 {V 6 - 6 {r ) - V|' 6 ], where the superscripts 
refer to the respective (3 values, and add this non- 
perturbative difference to the (3 = 6.6 aVs estimate 
above. In doing this we arrive at aVg P w 0.83 which 
we subtract from both NP and tbNLO data sets. Now 
at large R the tbNLO coupling will diverge linearly to- 
wards negative values as our V"J^ P estimate differs from 
Vg by as much as 0.35 a -1 . It is obvious when com- 
paring for instance the dimensionless NP couplings at 
R w 3.5 with R w 8.5 that both NP data sets scale 
nicely. We see however that at short distances the agree- 
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ment with the perturbative expectation at j3 = 6.0 is 
somewhat worse than that at (3 = 6.6 above. The gap 
between NP and tbNLO curves opens up in the same re- 
gion in lattice units, around R ss 4, which corresponds to 
a larger distance in physical units. Since the perturba- 
tive expectation decreases at large r only because of the 
difference between aVg F and aVg hNLO we should not ex- 
pect agreement at R > y5 anymore and interpret this 
behaviour as accidental. This view is confirmed by our 
j3 = 6.2 and (3 = 6.4 data sets: the improved agreement 
at small r goes along with the gap between perturbative 
and NP potentials opening up earlier. 



3. Comparison with rif — 2 

We compare the quenched potential with a potential 
obtained with sea quarks for the example of two flavours 
of Wilson fermions at (3 — 5.6 and k — 0.1565. These 
parameter values have been chosen to correspond to a 
quenched lattice coupling of (3 = 6.0: in the rif = 2 
case one finds [[39] r = 5.28(5)a while a quenched (3 = 
6.0 yields tq = 5.33 (3) a. The effective string tensions, 
determined by three parameter Cornell fits agree within 
errors too: era 2 = 0.0466(14) vs. era 2 = 0.0479(7), such 
that the linear term vanishes when subtracting the two 
potentials from each other. The k value corresponds to 
a quark mass ma ss 0.040 and we will use this value to 
obtain the perturbative expectations. 



1 
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Figs. |32] and [33] which is why we refrain from including 
LO, bLO, bNLO and tbLO curves. The respective NNLO 
results are not known at present, even at large R. 
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FIG. 38: The non-perturbative potential for Wilson fermions 
at P = 5.6, k = 0.1565, i.e. ma w 0.04, in comparison with 
NLO tadpole improved boosted and bare perturbation theory. 



In Fig. [38] we compare the NP rif — 2 potential with 
NLO and tbNLO perturbation theory. In the latter case 
we used the value Ajj^tq — 0.570(6), as obtained from 
the plaquette □ = 0.57073. The situation is very similar 
to that of the quenched potential at (3 = 6.0 depicted in 




FIG. 39: The difference between the n/ = 2 lattice poten- 
tial with Wilson fermions at (3 = 5.6, k = 0.1565 and the 
quenched potential at f3 = 6.0: non-perturbatively vs. per- 
turbatively. 

In a next step we calculate the difference between the 
rif = 2 and the n/ = potentials. This is depicted in 
Fig. [39], together with the respective differences calcu- 
lated in tbNLO and NLO. The NP linear confining term 
cancels in this combination, such that there is some rea- 
son to believe that this difference should be more per- 
turbative than the potentials themselves, at least as long 
as R < 12, the distance at which the QCD string will 
eventually "break" in the rif = 2 case. In the NP as well 
as perturbative data sets lattice artefacts are clearly vis- 
ible that can be ascribed to the effect of the sea quarks. 
In the case of NLO we have depicted the rif = expec- 
tation using the slightly smaller one loop matched value 
f3 = 5.969 rather than (3 — 6.0 for internal consistency. 

As R — > 00 all three curves approach constant values 
as the massive fermions start to decouple from the run- 
ning. Although the shape of all curves is similar, even the 
tbNLO prediction underestimates the NP data by almost 
30 %. This is not surprising as the vertical scale has been 
vastly inflated after subtracting two quantities of similar 
size from each other. For instance we can almost com- 
pletely eliminate the gap between the NP and tbNLO 

(2) 

results by increasing the AljL.r value (or decreasing the 

A^Lro value) by only 3 %. We find it quite encouraging 
that at least on a qualitative level we are able to resolve 
sea quark effects that only appear from 0(a 2 ) onwards 
in perturbation theory. Unsurprisingly, a quantitative 
understanding is not possible without knowledge of the 
difference to second non-trivial order, i.e. 0(a 3 ) in this 
case. 

We now turn to the difference in the running of the two 
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FIG. 40: The same as Fig. [B9| but with the difference 
of the respective self energy contributions subtracted, and 
multiplied by — r. We estimate aAV^ fa 0.062 and find 
aAVs LO « 0.042 and aAV| ,LO « 0.034. "3 tbNLO" refers 
to the tbNLO result, multiplied by three. 



respective Coulomb couplings, between rif — 2 and rif = 
and depict the quantity Ae(i?) = -RalV^ (Ra) - 
V"(°) (i?a) - AU S ] in Fig. |o| This combination further 
inflates scale and errors, in particular at large R. In the 
tbNLO case we can now identify deviations between the 
curve that constitutes the large R limit and the points 
calculated in lattice perturbation theory at short dis- 
tances. This bias, which also appeared in Fig. [l6] and 
which turns out t o be a n 0(ma) lattice artefact, has been 
discussed in Sec. Ill C above. Since at this lattice spac- 
ing, and within the scale of the figure, perturbative pre- 
dictions and NP data cannot be expected to agree with 
each other anymore (cf. Fig. [39]) we subtract different 
AVs values from the three data sets. While we compute 
aAV^ LO and aAV^ bNLO in perturbation theory, we fit 
aAU? p « 0.062 from the NP data depicted in Fig. ||. 

The NP data indicates an enhancement of the Coulomb 
coupling in the rif = 2 case at the given mass by al- 
most 0.03 over the quenched coupling at short distances. 
This enhancement then reduces to zero at large distances 
where the two data sets are matched. Contrary to this 
we find the bare NLO expectation to be rather flat as 
a function of R: the running of the coupling with R is 
faster when fermions are included, however, this effect is 



compensated for by a larger value of g 2 , in comparison 
with the quenched case. This larger bare coupling can 
be translated into a deflation of the distance scale which 
in turn slows down the running again. In the end one 
would expect an increase of e by 0.012 - 0.013 to NLO. 
The tbNLO prediction includes some slope but falls short 
of the NP data by a factor of about 3. After multiply- 
ing the points by this factor (open circles) the agreement 
with the NP Ae(R) is at best qualitative. We should not 



forget however that in the figure we have vastly inflated 
tiny differences. Unfortunately, rif — 2 lattice simula- 
tions at finer lattice spacings are not yet available such 
that we cannot test to what extent the quality of the 
parametrization improves as the continuum limit is ap- 
proached. 



V. SUMMARY 

We have calculated Wilson loops, the static potential 
and the static source self energy to 0(a 2 ) in lattice per- 
turbation theory with the Wilson gluonic action as well as 
with massive Wilson-Sheikholoslami-Wohlert and Kogut- 
Susskind fermions for arbitrary representation of the ex- 
ternal sources. Some numbers for "improved" gluonic 
actions are provided too. The results are useful for un- 
derstanding violations of rotational symmetry and the 
effect of including sea quarks on the lattice. Lattice re- 
sults can be related to other regularization schemes like 
the perturbatively defined MS scheme, for instance to 
allow for an extraction of the QCD coupling at high en- 
ergies from low energy hadronic phenomenology. 

We find that at presently employed quark masses lat- 
tice spacing effects on the short distance potential are not 
smaller in absolute terms for the two non-perturbatively 
0(a) improved fermionic actions, relative to the Wilson 
action. However, improvement results in a reduced mass 
dependence around ma = of quantities like the static 
source self energy, that is required for determinations of 
the b quark mass in an effective field theory framework, 
small Wilson loops and in the matching between lattice 
and MS schemes at finite a. 

We perturbatively determined (3 — 2N/g 2 shifts be- 
tween the pure gauge theory and QCD with massive sea 
quarks. The results appear to be reliable within 10 % for 
quark masses down to 20 MeV for Wilson type fermions 
(which translates into less than 1 % uncertainty in the 
value of the predicted matched 0) and within 25 % for 
"n / = 2" KS fermions within a window of lattice spacings 
that correspond to quenched f3 values 5.8 < (3 < 6.15. 
This method might turn out to be useful in predicting the 
impact of a shift in the sea quark mass onto correlation 
lengths like r$/a in future simulations. In particular we 
confirm that for SW fermions the lines of constant ro/a 
and m^a are significantly tilted with respect to those of 
constant (3 and ma, respectively, a fact that complicates 
continuum limit extrapolations. 

The success of the perturbative matching of the run- 
ning coupling at flavour thresholds via an intermedi- 
ate mass-dependent scheme, in our case the potential 
scheme, in lattice simulations (which allow us to compare 
with the correct non-perturbative result) means that the 
standard procedures employed in perturbative calcula- 
tions ^8) should be quite reliable even at the charm 
quark mass. 

We investigated so-called boosted perturbation theory 
methods (?], |^| and found them to work well in many 
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cases. However, such methods have to be digested with 
some caution: as boosted 0(a 2 ) perturbation theory is 
tuned to reduce higher order and in particular 0(a 3 ) cor- 
rections, differences between results calculated at these 
two orders are not necessarily indicative of the system- 
atic uncertainties involved. In the absence of even higher 
order calculations, the perturbative uncertainty in for 
instance predictions of the running coupling can be es- 
timated by varying both, the quantity from which the 
coupling is determined and the lattice spacing a. 

In analysing quenched lattice data on the static poten- 
tial and plaquette we obtain the value A^Lr = 0.61(9) 

which is in agreement with A^Lro = 0.60(5) obtained 
by use of finite size techniques by the ALPHA Collab- 
oration p35j. The error which is dominated by the un- 
certainty of higher order perturbative corrections can be 
systematically reduced by an 0(a 3 ) calculation of the 
static potential. Using Wilson, KS and SW n/ = 2 lat- 
tice data, the latter with Wilson as well as with Iwasaki 

(2) 

glue, we estimate Ai^L-ro = 0.69(15) in the continuum 
and chiral limit, where the error is again dominated by 
systematics. We find this value to depend strongly on 
both, the quark mass and the lattice spacing. Chiral and 
in particular continuum limit extrapolations are not fully 
under control but the use of different input quantities and 
lattice actions allows us to estimate the systematic un- 
certainties. 

Based on the n/ = and nf = 2 estimates we 

(31 

guess — 270(70) MeV. This value corresponds to 

a^L(m^) = 0.1133(59). In spite of the fact that more 
lattice data entered our analysis than has been used in 
any of the previous attempts to determine the n/ = 2 
running coupling, our error estimate is not smaller but 
larger than those that have been quoted by various lattice 
groups in the past @, §, §, g, g §, §§. The 

main reason 

is that the use of different quark actions and observables 
for the first time enabled us to realistically assess the 
uncertainties connected to higher order perturbative cor- 
rections and chiral and continuum extrapolations. 

We estimate that the error on the QCD coupling at 
the Z boson mass can be reduced to about 4 %, once a 
set of rif = 2 + 1 lattice simulations, similar to those that 
are already available for n/ = 2, becomes available. By 
improving the perturbative calculation to NNLO as well 
as by incorporating a wider range of lattice spacings into 
the numerical simulations, errors of 1 - 2 % appear to be 
realistic within the next few years. It might be possible 
to further improve on this quality by a dedicated lattice 
simulation of running quark masses and coupling p6| | . 

We found that the lattice potential at distances smaller 
than about 1 GeV -1 is quite well described by NLO 
boosted tadpole improved perturbation theory without 
any free parameter, as long as no attempt is made to sub- 
tract the power divergent static self energy. This obser- 
vation is somewhat in contrast to previous beliefs po| . |6^] 
and supports the view that perturbation theory for this 



quantity (and not only for the force) might be quite re- 
liable in the continuum too. This is so as long as renor- 
malons (in the MS scheme) and power divergences (in 
the lattice scheme) are properly dealt with when rear- 
ranging the perturbative series [^0[ At distances 
larger than 1 GeV -1 the linear confining term sets in 
rather rapidly and it is doubtful that this behaviour can 
be emulated completely by higher order perturbative cor- 
rections. At distances shorter than 1 GeV -1 the dif- 
ference between non-perturbative data and perturbative 
expectation is dominated by a quadratic term with a co- 
efficient of 0(1) in units of Ajj-g. The dependence of this 
coefficient on the way in which the perturbative series is 
organized, on the A^gTo value that is used and on the 
lattice spacing deserves further study. 

It turned out to be possible to resolve the difference 
in the logarithmic running of the coupling between two 
flavour QCD and the quenched approximation from non- 
perturbative data on the static potential. This difference 
is in rough qualitative agreement with the NLO expec- 
tation, however, about three times larger than expected. 
Unfortunately, at present we are limited to inverse lattice 
spacings of 2 GeV such that it is not possible to judge 
whether the discrepancy reduces as the scale is increased. 
As the magnitude of the effect should also increase as the 
quark mass is reduced, next generation lattice data at 
smaller quark masses should turn out helpful too. 
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APPENDIX A: PERTURBATIVE RELATIONS 
BETWEEN LATTICE AND CONTINUUM 
SCHEMES 

We review the known perturbative relations between 
the lattice /3-function and the modified minimal subtrac- 
tion (MS) scheme for various lattice actions in SU(N) 
gauge theory with nj fermion flavours. We also display 
continuum perturbation theory results for the (mass de- 
pendent) static potential in position space, which consti- 
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tute the limit against which our calculations will converge 
for r S> a. 



1. Conventions 



We define the QCD /3-function as, 
da 



P(a) 



din fi' 2 



-/3 a - /3iar - /3 2 a 



(Al) 



While in mass-independent schemes like lattice regular- 
ization 12 or minimal subtraction the first two coefficients, 
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(47T) 



(A2) 
a , (A3) 



are universal, /3 2 will in general depend on the renormal- 
ization scheme. In the MS scheme /3 2 reads K54], 



MS 
2 



2857 ,r3 


/1709 


N 3 - 
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V 54 
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187V/ 



iV" - 



187 



1 



36 AN 2 



7 



1 



(4tt) 3 



»/ 
(A4) 



The QCD A-parameter in a given scheme is defined as, 
1 



A = lim /i exp 

— *oo 



[Poa(n)] "* . (A5) 



We discuss the conversion between two different renor- 
malization schemes: let, 



a'{p) = a{n) + cia 2 (fi) + c 2 a 3 (^) + ■ 



(A6) 



From the definition of the /?- function, Eq. (Al), one in- 
fers, 

# = A), # = (A7) 

# - A - cyfa. + (c 2 - c 2 )/3 , (A8) 
A' = Ae Cl/(2M . (A9) 

If not only the scheme but also the matching scale n is 
shifted, one easily obtains the relation, 



«V) = a(jji) + cX(/i) + 4a (p) 



(A10) 



12 Strictly speaking, in lattice regularization, /3; only become scale 
independent in the continuum limit a — » 0, i.e. when scales 
fi <C a -1 and quark masses nii <C a -1 are considered. In gen- 
eral Pi will obtain additional non-universal contributions that 
are functions of ami and afi. While contributions of the latter 
type are genuinely non-perturbative in character, quark mass 
dependent contributions arise in perturbation theory too. 



with, 



ci - 2/?o In ( Z- ) , 
c 2 - (2ft + 4 Cl /3 ) In ( - 



+4/3 2 ln 2 (^ 



(All) 



(A12) 



c 2 - 2j8i In [ — | 



2 -c 2 
l c i- 



2. Conversion between the MS 1 and lattice schemes 

The MS coupling is related to the lattice coupling 
a L = g 2 /(4ir) via, 

a Ms( a 1 ) = a L + ha L + b 2 a 3 + ■ ■ ■ , (A13) 

where for the Wilson and Sheikholeslami-Wohlert (SW) 
quark actions one obtains M, |35|, |37], [7§, [73|] , 



6x 



2A 



+ kiN + K 1 (ma)n f , (A14) 



where K\(ma) = ^i(O) + AKi (ma), with the numerical 
values, 

fci = 2.135730074078457(2), (A15) 
Kx(0) = -0.08414443(6) 

+0.0634188788(l)c svy (A16) 
-0.3750240693(1)4^, 



for the Wilson parameter r = 1. 
of the SW coefficient, c$w = H 



jo) , (i) 2 , . . . 5 (o) 
c sw ' c swy ' ' lb c sw 
Susskind (KS) fermions Ki(0) reads fl37j |3j 

s (0) = -0.03298341916(1) 



The tree-level value 
3.3414(9)a L -I- • • • = 



1. In the case of Kogut- 
(A17) 



Note that in our notation rif KS quark flavours corre- 
spond to rif (not 4n/) flavours of Dirac quarks in the 
continuum limit. Therefore, n/ has to be a multiple of 
four. The value of k\ applies to the Wilson gluonic action 
only. However, the contribution Ki(ma) (whose quark 
mass dependence has been calculated for the first time 
in this publication and is displayed in Tab. |x|) remains 
the same for Symanzik type improvements of the Wilson 
gluonic action. 

At vanishing quark mass or lattice spacing, m <C a -1 
we obtain the massless limit, i.e. AATi(O) = 0. The 
mass dependent term results in a change of the QCD 
/3-function at finite a: 



At = Af s - 



a dAKi(m,ia) 



^ 2 



da 



(A18) 



From this we can infer, \m\ ma ^ dAKi (ma) /d In a = 

0: the dependence is 



lim ma _ + o dAKi (ma) / d In (am) 
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regular in ma near the origin. Since the /3-function is no 
on-shell quantity, we cannot exclude terms proportional 
to ma, even for improved quark actions. By demanding 
that in the limit ma — > oo the quenched /3-function is 
re-obtained, we find, 



lim 



dAKi(x) 
dlnx 



1 

3^' 



(A19) 



i.e. — AK\ will grow logarithmically for large ma. This is 
expected since ma — » oo means that the coupling on the 
left hand side of Eq. ( fA13| ) runs with nt active flavours 
while the running on the lattice side of this equation cor- 
responds to that of pure Yang-Mills theory. The differ- 
ence has to be compensated for by AKi(ma). We will 
determine the integration constant in the region ma — > oo 
in Eq. ( B19| ) below from matching the un-quenched and 
quenched static potentials in the infra red. 

The coefficient b 2 is at present only known for massless 
Wilson-SW quarks and Wilson gluonic action: 



3tt 2 
87V 2 



+ k 2 + k 3 N 2 + 



K2 
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(A20) 



with 1 35, M 



k 2 
h 
K 2 



From Eqs. 



-2.8626215972(6), 

1.24911585(3), 

0.1890(2) 

-0.02492(5)c SH / - 0.83585(3)c| w 
-0.79942(5)4^ - 0.012947(2)4^, 
-0.1579(3) 

-0.00540(6)c 5W + 0.7684(l)4iy 
+0.033842(5)4w + 0.006920(2)4^. 

and (|A9|) we infer, 
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This results for instance in the numerical values, 
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MS 

for N = 3 colours 
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(A27) 
(A28) 



Conversion between the MS scheme and 
schemes based on the static potential 



where q = |q| and Cf — (N 2 — 1)/(2N). The momentum 
space potential is related to that in position space via a 
Fourier transformation, 



V(q) = 



d 3 x 



: [V(r)-V s ]e^, (A30) 



where r = |x|. Vs denotes the self energy term, that 
vanishes in dimensional regularization but will in general 
be present, and the position space potential, 



V(r) = -C F aR{r l) +V S , (A31) 



defines the coupling an. By Fourier transformation one 
obtains the relation between the two schemes, 



OLR 



(M) = oty (e-» (l + ^L 2 ) 



(A32) 



which holds for fermion masses m, <C /i. 
0.57721566 ... is the Euler constant. 

We write, 



av(jJ>) = °<-ms{v) + aiaj^(n) + a 2 a 3 H . (A33) 



The coefficients a\ and a 2 are known 13 . For nf massless 
quark flavours they read (33J : 
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(A35) 



We define the coupling a v through the static QCD ^ denotes the value of the (-function C(3) = 

potential in momentum space, 1.20205690 . . .. The corresponding result for massive 

1 \ quark flavours with masses mi has been obtained in 

V(q) = -4irC F ^f-, (A29) Ref. @: 
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13 In fact the leading log contribution to CI3 ^3 and, more recently, 
the complete NNLO result |75| have been calculated too. 
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with the numerical values Il4l| 



Co = 5.19(3), 

d = -0.571(34), 

C 2 = 0.221(15), 

C 3 = 1.33(12). 



(A38) 
(A39) 
(A40) 
(A41) 



A c hange in the argument on the left hand side of 
Eq. ( A33| ) , ay(/i) 1— > av(fJt'), results in the substitutions, 



oiOi'.W) = ai (imi^jj -2/3 m(^j , (A42) 
Oa0*',{mi}) = a 2 fjm^jj +4/3 2 ln 2 fej (A43) 



2/3 1 +4/3 ai [\ m ijj 



in analogy to Eqs. (All) - (A12) 



The masses m, in Eqs. ( |A3q ) - (|A37|) denote pole 
masses, 



Cf 

47T 



4 - 3 In ■ 



"W/*)- (A44) 



Re-expressing Eq. ( A37 ) in terms of MS quark masses 
rrii at a scale fj,, therefore, results in the replacement, 

2 

n ^ 4 3 In m * 

°P(W) = aaCW) + ^ E " / > (A45) 

i=l + Com? 



while Eq. ( A36 ) remains unaffected. For Wilson fermions 
one obtains~i|76| , 

TO Ms( a_1 ) = ( 1 + 16 - 95241 ^- a ) m o- ( A46 ) 



Hence, in terms of lattice quark masses TOj(a), a 2 reads: 

C f ^4 16.95241 - 3 In mV 



a 2 ({™*» = a 2 ({m l }) + Y^a" E 



£ 1 + (C m 2 a 2 ) 



We now go to position space and write, 



(A47) 



a R (n) = «ms(m) + afaj^ip) + a£a 3 + ■■■. (A48) 



In following Eqs^jARj - ( |A13| ) we obtain from 
Eqs. ( |A3l ) and 



af = 01 + 27/%, 

= a 2 + 47/? ai + 27ft + ( 47 



(A49) 



2 1 y )/?o 2 (A50) 



The one loop result for massive quark flavours reads 
af({m i }) = af + i-^Ein^ (: "'" 



with 



Ein(x) = 7 + ln(ar) + J5 X (or) = - V" 



{- x y 



(A51) 



(A52) 



For the corresponding two-loop mass dependence of the 
position space potential we refer to Ref. p4| . The expo- 
nential integral E\(x) is defined by, 



dt- 



f x 1 - e"* 
Ein(x) = / dt . 



(A53) 



(A54) 
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From Eqs. flA6j), ( |A9| ), ( |A33| ) and ( |A34| ) one obtains 
the numerical values, 



A<?>« 1.5995 A^. 



A 



(2) 



while Eq. (A32) results in 



(«/) 



>A 



(«/) 



1.5213 A^, 



1.7811 A 



("/) 
v ■ 



(A55) 



(A56) 



Note that for massive quarks the coefficients of the V- 
and i?-scheme /3-functions become explicitly scale de- 
pendent. However, this does not affect the respective 
QCD A-p arameters, that are defined in the limit /i — *■ oo 
[Eq. ©]. 



We obtain: 



dan 
d In /i 2 



day 
din ji 2 



ft 



(«/) 



(»/) 



n f 

-E 



n f . 



()7T 



1 — exp — 




(A57) 
(A58) 



In the situation, mi = m 2 = ■ ■ ■ = m n ._i = 0, m nj > 0, 
the above coefficient functions smoothly interpolate be- 
tween /3g™ f 1 ' (for /j,/m n , — ► 0) and /3q (for fi/m nf — > 
oo) as one would expect. 



4. The perturbative potential at large i?. 

We denote a dimensionless distance measured in lattice 
units by R while r = Ra refers to the distance in phys- 
ical units. At separations R 1, i.e. r 3> a, rotational 
invariance should be restored and the static potential ob- 
tained in lattice perturbation theory should coincide with 
the one computed above by use of dimensional regulariza- 
tion. In the large R limit our results should therefore ap- 
proach the continuum expressions given in Appendix A 3 
above. 

We define, 



aV(Ra) = Vl (R)a L + v 2 (R)a 2 L + ■ ■ ■ , (A59) 
aVs = vi(oo)a L + v 2 (oo)a 2 L H , (A60) 

where the self-energy Vs = 2<5m sta t is twice the lattice 
pole mass of a static colour source. The coefficient func- 
tions i>i(R) have been calculated in the present article. 
At large separations we expect Uj(R) only to depend on 
the modulus R = |R| and can write, 



v c ,i(R) = Vi(R) - Uf(oo). 



(A61) 



Eq. (A31) now implies, 
a [V(Ra) - Vs] — ► 
It follows, 



an [{Rar 1 ] 
-Lf —p, (it —> oo) 



R 



C F 



(A62) 



«i(R)-»i(oo) — »« c ,i(i2) = -- f (R^oo), (A63) 
where 



for the Wilson gluonic action. 12(00) is also known. Here 
we state the result for massless Wilson-SW quark flavours 
and N = 3: 

i 2 (oo) = C F x {16.714(4) + [-0.42333(6) (A65) 
+0.0516(2)c 51 y - 0.5870(2)4^71/}. 

In the case of massless KS quarks the numerical value 
of the rif coefficient is — Cf x 0.36846(6). Results on 
the mass depe nden ce of the fermionic coefficients can be 
found in Tabs. XII and XII] . The rif = result has first 
been derived in Rcfs. ]13j, |15|. The precision has subse- 
quently been increased by Martinelli and Sachrajda 
We do, however, slightly disagree with this reference and 
obtain 5m s tat,2 = Cfv 2 (oo)/2 = 11.143(3), rather than 
11.152 ||. The latter reference includes the Wilson and 
SW contributions for massless fermions as well, which 
we were able to reproduce with increased precision. The 
mass dependence of 12(00) as well as the KS result are 
new. The rif = value for 113(00) has recently been ob- 
tained by means of stochastic perturbation theory [plj : 



i 3 (oo) = C F x 129(2). 



(A66) 



This result agrees with the one of Ref. p2] 13(00) = Cf x 
130(1), obtained from a fit to high (3 Monte Carlo data. 
From Eqs. flA59| ) - ( |A63| ) we obtain, 



a R [(Ra) x ] = a L + ei(.R)a£ + e 2 (R)a 3 L 
with 

ei{R) = T^r- = -— i c , l+ i i? • 

v cA (R) C F 



(A67) 



(A68) 



11(00) = C F x 3.1759115. 



(A64) 



In converti ng th e scheme and scale by use of Eqs. ( A4S| ). 
( |Al3j ) and ( |A10| ) - ( |A12| ), we obtain the coefficients ejjR) 
and e 2 (R): 
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a R [(Ra)- 1 ] = aj^KRa)- 1 } + of ({Rm^a^iRa)- 1 ] + a^({R mi })a 3 

= a L + [bx + af ({Rrrn}) + 2/3 lni?] a 2 + {b 2 + af ({Em*}) + 2b 1 a*({Rm l }) 

+ [2/3i +4(6i +a?{{Rm l }))Po}lnR + 4(3o\n R}a 3 . (A69) 



Therefore, 

v c , 2 (R) = 

M^) = 
This amounts to 



R 
C F 



[61 + af ({Rrrn}) + 2/3 In JZ] , (A70) 



{fc + af ({i?TO,}) + 26iaf ({i?™,}) + [2/?i + 4(&i + af ({i?TO,}))/3o] In R + 4/3 2 In 2 i?} . (A71) 



C 



, 2 (R)=-^-\A +AilnR + Y / 



AKi(rriia) + — Ein(vCoTOiai?) 

07T 



(A72) 



with 141516 the numerical values for n/ flavours of Wilson- 
SW fermions in S77(3) QCD, 

A 7.7164259753 (A73) 
- [0.233808 - 0.063419(l)c SV y + 0.37524(l)c| w > / , 
Ai ps 1.75070437401- 0.10610329540 n f . (A74) 

For KS fermions the fermionic contribution to Aq differs: 

Aq S ps 7.7164259753 - 0.1826473162 n f . (A75) 

The function Ein(x) is defined in Eq. ( A.52Q and can be 



14 We have factorized out the Casimir constant Cp such that the 
above results also apply to the potential between charges in rep- 
resentation D, replacing Cf by Co. For TV ^ 3 the constant Co 
has not yet been calculated. 

15 Rather than comparing the perturbative lattice po tentia l at large 
distances with known results, one can also use Eq. ( A70 ) to deter- 

hl(H) - v i+1 (oo) 
a no way 



mine the parameters bi from the calculated Uj. 
at large R. While in precision this procedure can 
compete with the values displayed in Eqs. (A.16) or (A17), cal- 
culating big Wilson loops and extracting the attice potential by 
means of stochastic perturbation theory [3ll 77] might turn out 
to be a feasible alternative to diagrammatic techniques for more 
involved quark and gluonic actions. 

v c i are subject to (am)" and (a/r) 1 ' lattice corrections, where 
v = 1 in the case of Wilson fermions and v = 2 for the KS and 
SW actions. While the r dependent corrections vanish at large 
r where we match to the continuum scheme, the am corrections 
will in general be present in this limit. For < ma < 00 AKi can 
be defined by matching the lattice scheme to a mass-dependent 
scheme like the R-scheme. We do this by imposing that AKi con- 
tains all am corrections to v Ci 2- One can use different quantities 
to define AKi, for instance the potential in momentum space. 
Since in general the (ma) u corrections for different observables 
will differ, AKi can only be fixed up to an (ma) v ambiguity. 
This is a reflection of the loss of universality of the /3-function, 
due to the mass dependence of /3o (and Pi). In the continuum 
limit, a — > universality will be restored. 



approximated for small arguments by, 
— Ein ^y^CoTOjai?^ = ^J~CQrriia + [(m^ai?) 2 ] 



(A76) 



to leading order the coefficient v Ct 2(R) is parallel shifted 
in proportion to ma at short distances. 

At infra red scales, Ra 3> m^ 1 , E\(x) oc e~ x ' /x, such 
that the massive flavours decouple from the logarithmic 
running of the potential. For nf massive degenerate 
flavours with masses to,- = to, wc find: 



„(«/) 



(R) 



u c,2 



(R) 



{*<»»> +s 



(A77) 



6 



For massless fermions we can parametrize the next or- 
der t erm i n the expansion of the lattice potential v c 3 by 
[Eq. ([A7j])], 



V c ,3 



Cf 
R 



Bq + Bx \nR + B 2 In 2 R) 



(A78) 



The fermionic part of 62, which is required to calculate 
the contribution to Bq that is proportional to nj, is only 
known for massless Wilson-SW fermions. Restricting 
ourselves to massless Wilson quarks and SU(3) we ob- 
tain the numerical values, 



B 



Bo 



73.8166085803 (A79) 
-3.667(1) n f + 0.0530952 n), 
28.3102065046 (A80) 
-2. 295710a/ + 0.0496157 n), 
3.06496580518 (A81) 
-0.371511006771/ + 0.01125790929 n^. 



B 2 is known exactly while the accuracy of the term within 
B that is proportional ton/ is limited by the precision 



of the constants K 2 and K 3 [Eqs. (A23) - (A24)] 



42 



APPENDIX B: THE PERTURB ATIVE /3-SHIFT 



1. Matching quenched and un-quenched 



In Eq. ( A77 ) we have seen that at distances r ^> m~ l 
the running of the coupling is not affected by the pres- 
ence of sea quarks anymore: at large distances, at least in 
perturbation theory, the effect of massive quarks can be 
integrated out into a shift of the coupling constant of the 
quenched theory. The matching can be done within an in- 
termediate mass-dependent scheme, such as the V- or R- 
schemes. For simplicity we assume ri/ mass-degenerate 
flavours. The discussion below can easily be generalized 
to the non-degenerate case. For Ra 3> mT x we write, 



a[V(Ra) - V s ] 
v e ,i(R) 



^+eft\R)al 



„(°) _l .(o) 



r\R)ai 



(Bl) 



where ef 1 ^ (R) are the expansion coefficients of Eq. ( A67 ) 
for rif massive flavours with lattice quark mass m = 
— K~ 1 )/(2a). The matched coupling that results 
in the same physics at R (ma) -1 is, 



,f = Z(m)a^ 



with 



Z(m) = 1 + zi(m)a^ s ^ + Z2{m)a 2 L + ■ ■ • , 



z\{m) 



lim 

R^ao 



e^ f \R)-e^{R) 



(B2) 



(B3) 



(B4) 



z 2 (m) = lim \e^ f) (R) 



(R) 



-2e 



?\R)[e^(R)-e^(R)}}. (B5) 



In fact it turns out to be easier to determine Z(m) in 
momentum space in the limit q — > 0. Since one and 
the same physical potential is matched, momentum and 
position space results have to agree [modulo the (am) 1 ' 
ambiguity within AK{\ and we explicitly checked this in 
our calculation. 

Using the MS scheme in momentum space we demand 
for q — > 0: 



4 lf \l) = A ) ( a - 1 )+4" /) (9,{m}) a 2 + ... 



MS 
MS 



(a- 1 )+af\q)a 2 + 



(B6) 



where the ai(q,{m}) are defined in Eqs. ( A42j ) - ( A43| ). 
This amounts to, 



(0) 

t— — 

MS 



(a-')=W(m)a^(a-') 



(B7) 



with 



W(m) = l + Wl (m)a^ (a" 1 ) + w 2 (m)a^ (a" 1 ) +• ■ ■ 



(B8) 



wi(m) 



lim 

q^0 



a ( ^\q,{m})-a^(q) 



(0), 



(B9) 



a , - a , ■ + In (^fc^ma^j , 



,(»/) 
n 



AO) _l «/ 



w 2 (m) = lim \a { 2 f \q > {m})-a 2 0) (q) (BIO) 

q^0 I 

— 2(4° (q)wi(m) 



By use of Eq. ( Al3| ) one can then relate the above MS 
scheme coefficients to the lattice coefficients of Eq. 



z 1=Wl + b ( ™ f) - bf } =wi + Kw f , 



Z-2 



(Bll) 



(0) 



W2 + b ^- b W +2 (b^-bW) ( Wl -b[ 
w 2 + + K 3 N + 2K lWl + K 2 nA n f . (B12) 



2. Results 

In following the above considerations we obtain the one 
loop result, 



Z\(m) = 



5 \ / 

Ki(ma) - — — + — In ( V C ma 
18n Sir \ 



the infra red behaviour of the theory with massive 
fcrmions is the same as that in the absence of massive 
flavours, at a gauge coupling that decreases with the 
quark mass: 



a 



(0) 



a 

7) 

3tt 



[In(Dma) + 3TrAK 1 (ma)} a 2 , 



D 



^Coexp[3-KKx(0) - 5/6]. 



The constant D has the numerical values, 



D = 0.448(2), 



D 



sw 

^KS 



0.0238(1), 
D KS = 0.726(2), 



(B14) 



(B15) 



(B16) 
(B17) 
(B18) 



for Wilson, SW and KS fermions, respectively. From 
the values of D we see that the coupling shift is much 
more pronounced for SW fermions than it is for Wilson 
or KS flavours. From Eq. (Al£) we infer that for large 
ma Aifi (ma) = -ln(D'ma)Jl3ir). Note that Eq. ( |bT^ ) 
fixes the coefficient D 1 = D and therefore the leading 
order behaviour of AK i (ma) : 



AKx(ma) 



1 

3^ 



ln(Dma) 



(B19) 



0.00105(30) - ifi(0) 



3tt 



In ma. 



In the above limit as well as in the limit am = 
AK\(ma) is uniquely determined. This is different in 
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the intermediate ma range where AKi (ma) will depend 
on the quantity used in the matching and is therefore 
subject to 0[(amY] ambiguities. This in turn will also 
affect the predicted /3-shift. 

The coefficient of the a 2 term diverges for m -C a^ 1 : 
for small values of m, the physics at distance scales 
r 3> to -1 , at which the matching to the quenched the- 
ory has to be performed, will be dominated by non- 



perturbative effects. Eventually, at very small quark 
masses the quenched and un-quenched theories cannot 
be matched to each other anymore at any distance scale 
with relevance to hadronic physics and the un-quenched 
theory becomes genuinely different from the quenched 
one. 

Finally we state the result for 2:2 (m): 



z 2 (m) = 



62 56 \ N 
" l T + T C3 j T~ 



38, (m 2 a 2 
— In — ==■ 
3 



1 



16tt 2 



100 4 l2 , _ 

V - In 2 (C r 

81 9 V 



19 C 3 -2C*i 



'16.95241 - 3 In to 2 a 2 



+ 



16tt 2 




(B20) 



where m denotes the lattice mass (k^ 1 — n~ x )/(2a), 
rather than the pole mass (hence the correction term con- 
taining the numerical factor 16.95241). Unfortunately, 
we cannot yet apply this result to lattice data. For even 

17 In the case of Wilson fermions one might be tempted to argue 
that 0(ma) terms are lattice artefacts anyway. However, in the 
determination of the /3-shift it is exactly this ma dependence 
of ro/a that we are looking for and this indeed depends on the 



in the case of Wilson-SW fermions where ^(0) and 
K 3 (0) are known, the quark mass dependence of these 
coefficients still has to be computed 17 . 

lattice action. Hence the variation of K2 and Kz with ma cannot 
be neglected. 
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